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We study theoretically the energy and spatially resolved local density of states (LDoS) in graphene 
at high perpendicular magnetic field. For this purpose, we extend from the Schrodinger to the Dirac 
case a semicoherent-state Green's-function formalism, devised to obtain in a quantitative way the 
lifting of the Landau-level degeneracy in the presence of smooth confinement and smooth disor- 
dered potentials. Our general technique, which rigorously describes quantum-mechanical motion 
in a magnetic field beyond the semi-classical guiding center picture of vanishing magnetic length 
(both for the ordinary two-dimensional electron gas and graphene), is connected to the deformation 
(Weyl) quantization theory in phase space developed in mathematical physics. For generic quadratic 
potentials of either scalar (i.e., electrostatic) or mass (i.e., associated with coupling to the substrate) 
types, we exactly solve the regime of large magnetic field (yet at finite magnetic length - formally, 
this amounts to considering an infinite Fermi velocity) where Landau-level mixing becomes negli- 
gible. Hence, we obtain a closed-form expression for the graphene Green's function in this regime, 
providing analytically the discrete energy spectra for both cases of scalar and mass parabolic con- 
finement. Furthermore, the coherent-state representation is shown to display a hierarchy of local 
energy scales ordered by powers of the magnetic length and successive spatial derivatives of the 
local potential, which allows one to devise controlled approximation schemes at finite temperature 
for arbitrary and possibly disordered potential landscapes. As an application, we derive general 
analytical non-perturbative expressions for the LDoS, which may serve as a good starting point for 
interpreting experimental studies. For instance, we are able to account for many puzzling features 
of the LDoS recently observed by high magnetic field scanning tunneling spectroscopy experiments 
on graphene, such as a roughly -^/m-increase in the mth Landau-level linewidth in the LDoS peaks 
at low temperatures, together with a flattening of the spatial variations in the Landau-level effective 
energies at increasing m. 

PACS numbers: 71.70.Di,73.22.Pr,73.43.Cd,03.65.Sq 



I. INTRODUCTION 

A. Quantum-Hall effect in graphene 

The observation of an anomalous quantization of the 
Hall resistance in graphene at high magnetic fields^— 
related to the massless, relativistic-like spectrum of low- 
energy electrons on the two-dimensional honeycomb lat- 
tice, has triggered much excitation in recent years, see 
Ref. |H for a review. Indeed, the experimentally mea- 
sured Hall resistance follows the Landau-level structure 
expected for massless Dirac electrons;^ E m = ±^fmMl c 
in the clean case, with m a positive integer and il c = 
V%vf/Ib the graphene characteristic frequency given in 
terms of the Fermi velocity vf and of the magnetic length 
Ib = \Jhcl\e\B (here e = — |e| is the electron charge, c 
the speed of light, and B the magnetic field strength). 
The \[~B dependence of the characteristic frequency fi c 
in graphene, to be contrasted with the linear dependence 
of the cyclotron frequency w c = \e\B/(m*c) of more stan- 
dard two-dimensional electron gases (2DEGs) based on 



semiconducting heterostructures (in this case, m* is the 
electronic effective mass) described by Schrodinger equa- 
tion, constitutes one of the main signatures used so far 
in experiments to exhibit the relativistic-like character of 
the massless charge carriers. 

Also quite remarkable is that graphene displays a sur- 
face opened to the outside world, providing a direct win- 
dow to its electronic excitations. This is a clear ex- 
perimental advantage of graphene compared to 2DEGs 
based on semiconducting heterostructures, where the 
2DEG is buried deep inside the structure (typically 100 
nm or more). Graphene thus offers the opportunity to 
obtain precise insights into local physical properties of 
quantum-Hall systems, such as the local density of states 
(LDoS) via scanning tunneling spectroscopy (STS) mea- 
surements. In contrast, such local probes experiments 
have very poor spatial resolution in ordinary heterostruc- 
tures, although some progress has been made recently, 
see Ref. 0- This technical advantage will be certainly 
important in the future to elucidate the relation between 
microscopic inhomogeneitics induced by various disorder 
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types and macroscopic transport properties of large sam- 
ples. Various open questions in this respect are the na- 
ture of the universal plateau to plateau quantum phase 
transition^— or on a more quantitative level the precise 
formation of wide Hall plateaus. To pursue this goal, 
STS is one of the interesting available experimental tech- 
niques, and first experiments in graphene at high mag- 
netic field have been performed recently] 11 ' 12 Since this 
spectroscopic method gives direct information on the lo- 
cal electronic states, a better understanding of the LDoS, 
specific to the case of graphene at high magnetic fields 
and in arbitrary potential landscapes (without proceed- 
ing to disorder averaging), needs to be achieved. This is 
the main aim of the present paper. A second important 
aspect of our work is to obtain analytical solutions for 
a large class of parabolic confinement models, and as a 
motivation we now discuss the different types of poten- 
tials that can be involved in the two-dimensional Dirac 
Hamiltonian. 



B. Disorder types for graphene 

Because of the multicomponent structure of the wave 
function for graphene, several types of disorder can oc- 
cur, which we introduce here. The quasiparticle disper- 
sion for graphene has two Dirac cones (two "valleys" ) at 
low energies. For a given valley, the Hamiltonian in the 
presence of a perpendicular magnetic field has a matrix 
structure and is written as 

H Q = v f <t • f[, (1) 

where vf is the Fermi velocity, <x is a vector whose com- 
ponents are the Pauli matrices a x and a y in the pseu- 
dospin space, and the momentum operator is 

ri = -ihV r - -AM. (2) 

c 

The vector potential A is related to the uniform trans- 
verse magnetic field B via the relation V x A = B = Bz. 
For convenience, we will omit both physical spin and 
valley indices, thus assuming that the two valleys of 
graphene remain completely decoupled from each other 
and can be studied separately.— 

Quite generally, potential terms appear as either a ran- 
dom scalar potential, a random Dirac mass or a random 
vector potential^ The Hamiltonian in presence of these 
potentials is given by 

H = H + V(r), (3) 

where the function V(r) takes the general form 

V(r)= a P V P^ W 

p=s 1 x 1 y 1 z 

with a s the identity matrix in the pseudospin space, as- 
sociated to the scalar potential term V^(r). This contri- 
bution may have many different physical origins: electro- 
static confinement potential, impurity random potential, 



and/or Hartrcc potential resulting from the mean-field 
mutual Coulomb interaction between the electrons. The 
diagonal but antisymmetric term V z (r), associated to the 
o z Pauli matrix, describes the so-called random mass po- 
tential. This contribution might be introduced by the 
underlying substrate in single-layer graphene, while in 
bilayer graphene, such a term can be produced in a con- 
trollable way by introducing different electrostatic poten- 
tials in the two layers*^ The off-diagonal contributions 
coming as V(r) = [V x (r) , V y (r)] can be associated with 
a random vector potential, coming from the spatial dis- 
tortion of the graphene sheet in the third dimension by 
ripples4iii In what follows, all three possible types of 
disorder will be considered within the high magnetic field 
regime. 



C. Existing theoretical results for graphene in 
various potential types 

Let us first discuss various toy models of potentials 
(in a magnetic field) that were studied in the recent 
graphene literature. Quite generally, within the Dirac 
equation fewer models can be solved exactly than within 
its non-relativistic counterpart. For instance, the classic 
one-dimensional parabolic confinement model, as well as 
the circular parabolic confinement model, are seemingly 
not analytically tractable. For the 2DEG, the former is 
the well-known model to introduce the edge states and 
explain the quantized conductance in Hall bars. The lat- 
ter is the basic model for quantum dots and leads under 
magnetic field to the Fock-Darwin states with discrete 
energies. Thus, only much simpler models can be solved 
analytically for graphene, such as the uniform electric 
field . 15 ' 16 Progress can be achieved for circular hard-wall 
confinement with either scalar— or mass^ potentials, but 
only a solution in terms of special functions is then pos- 
sible. For parabolic and more complex potentials, fully 
numerical methods have to be used, e.g., see Ref. h» We 
will show in this paper that the limit of negligible Lan- 
dau level mixing allows one to solve analytically a large 
class of parabolic models, providing new insights in the 
high magnetic field regime. 

Coming to the more complex question of disorder, even 
less is actually known. Recent work devoted to the 
quantum-Hall effect in graphene has proposed to take 
into account disorder phenomenologically in the expres- 
sion of Green's function by adding a constant imaginary 
part iT in the self- energy— £H but Hall quantization ob- 
tains only in the limit where the energy rate r — > 0. The 
LDoS in the vicinity of a single pointlike impurity and in 
the presence of a strong magnetic field has been studied 
recently . 21 ' 22 Various types of disorders were also consid- 
ered in Refs. [23| and within the self-consistent Born 
approximation. While this method may be justified for 
short-range scatterers, it turns oulj 2 ^ to be inappropriate 
for a smooth potential in high magnetic fields. Because a 
quasi-local picture takes place in the high magnetic field 
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regime^ our calculation will be able to provide accurate 
expression for the LDoS in smooth arbitrary potentials. 



D. High magnetic field regime 

The strategy to follow is best explained by starting to 
discuss the specific nature of disorder for 2DEGs at high 
magnetic field. For very clean heterostructures, the dis- 
ordered potential seen by the electrons is mostly smooth 
on large length scales (several tens of nanometers), as 
the majority of impurities sit far away from the 2DEG. In 
contrast to the low magnetic field regime, where the elec- 
trons explore ergodically macroscopic regions of the sam- 
ple, the high field regime is characterized by cyclotron 
motion close to equipotcntial lines of potential landscape 
V(r) with a narrow transverse spread proportional to the 
magnetic length (which is smaller than 10 nm at several 
tesla) . The disorder landscape felt by the electronic wave 
functions is therefore very smooth in that situation. We 
note that in graphene additional sharper potential vari- 
ations (such as atomic vacancies of the carbon layer, or 
local imperfections from the nearby substrate) can oc- 
cur, although these tend to be detrimental to quantum- 
Hall physics by increasing the mixing of Landau levels. 
The coupling to the substrate can however be removed 
by suspending graphene flakes or with a decoupled layer 
in epitaxial graphene^ resulting in very high mobility 
samples. In fact, for both non-relativistic 2DEGs and 
graphene, the essence of the quantum Hall effect lies al- 
ready by considering smooth potential variations only, 
which is the case to be followed from now on. 

Theoretically, this smooth disorder regime was shown 
to be problematic at high magnetic field for standard 
quantum-mechanical methods based on perturbative ex- 
pansions in potential strength.-^ In that case, the high 
magnetic field regime is the correct starting point, and 
is characterized by two different dimensionless small pa- 
rameters: (i) Ib/(, associated to the transverse spread 
of the wave function along the classical guiding cen- 
ter R, with Ib the magnetic length and £ the typical 
length scale related to local variations in the potential; 
(ii) Ib\ 'VV\/fiiJ c ~ IbSV/^HuJc associated to Landau level 
mixing by local gradients | VV| of the potential, introduc- 
ing SV the typical amplitude variations in the potential 
on the scale £, and the cyclotron frequency lo c in the 
2DEG case. 

Clearly, quantum mechanics calls for non-zero Ib/S,, 
otherwise the so-called semiclassical guiding center pic- 
ture at Ib = emerges, giving at best a qualitative 
picture, and missing important quantum effects such as 
level quantization, tunneling, or interferences effects due 
to the potential energy V(R). The second parameter 
Ib\VV\/?iuj c controls the degree of Landau level mix- 
ing, so that Landau levels strictly decouple at infinite 
lo c . Most previous works have considered either limits 
separately (either Ib — > or w c — > oo), and the necessary 
formalism to incorporate both non-zero Ib and finite w c 



was developed for the standard 2DEG by the authors in 
Refs. I2a - l28l which will be extended in the present paper 
to the case of graphene. This mathematical construc- 
tion shows that a local picture of the high magnetic field 
physics emerges in terms of semicoherent-state Green's 
function, with a hierarchy of local energy scales^ 6 - ordered 
by powers of the magnetic length and successive spatial 
derivatives of the confinement or disordered potential. 



In the simplified, yet fully quantum limit of infinite 
cyclotron frequency and non-zero Ib , initial progress was 
made by other authors in Refs. [29| and [13 for the 2DEG 
case, where it was shown that Schrddinger equation ac- 
quires a uni-dimcnsional character, offering an analysis 
for toy models of confinement or tunneling in the low- 
est Landau level. The general structure of this limit was 
clarified in further developments in the Green's-function 
formalism , 26 ! 31 and this will be also examined in detail 
for graphene in the present paper. Our methodology is 
based on the exclusive use of Green's functions, not wave 
functions, for the simple reason that we project the quan- 
tum dynamics onto a semi-coherent representation with 
nonorthogonal states, forcing us de facto to give up the 
wave-functions picture. Noticeably, because the over- 
complete character of the chosen representation allows 
one to get rid of the Hilbert-space formulation inherent 
to the traditional operator formulation of quantum me- 
chanics, a unification of closed and open systems quan- 
tum mechanics is made possible here, i.e., one can get and 
treat quantization and lifetime effects on an equal foot- 
ing. An important application is the possibility to write 
down in the 2DEG case a unique Green's-function ex- 
pression which holds for all cases of quadratic potentials. 
This derivation has clearly proved that the appearance of 
lifetimes (expressing the presence of decaying states, i.e., 
an intrinsic time asymmetry) has for physical origin the 
instability of the dynamics occurr ing at saddle points of 
the potential landscape; see Ref. \2Q for a thorough dis- 
cussion of this point. 



In the graphene case, our calculation at large char- 
acteristic frequency fi c (or cquivalently at large Fermi 
velocity) brings important information, because, in con- 
trast to standard 2DEGs, even simple models of parabolic 
confinement for graphene do not possess an analytic so- 
lution at finite fl c . However, we will show that the limit 
O c — > +oo is exactly solvable for most quadratic poten- 
tials, allowing us to extract the explicit discrete energy 
spectrum in case of several parabolic confinement mod- 
els, and also, in principle, the transmission coefficients 
in case of tunneling near saddle points. Going beyond 
these toy models, our general formalism also allows us 
to calculate in a controlled way the LDoS in an arbi- 
trary and possibly disordered potential landscape. Our 
results will be discussed with respect to recent experi- 
mental findings ! 11 ! 12 
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E. Structure of the paper and summary of results 

First, in Sec. HU we shall investigate the free Dirac 
Hamiltonian in a transverse magnetic field, and introduce 
the graphene vortex states, which are the building blocks 
of the whole theory developed here. These states form 
an overcomplete family of semicoherent states, strongly 
localized around arbitrary guiding center positions R, 
and encode the cyclotron motion quantum mechanically. 

In subsequent Sec. IIII1 we introduce the Green's func- 
tion for graphene vortex states and derive its general 
equation of motion, Eq. (|50j) , including Landau-level mix- 
ing processes. The general connection to the real-space 
Green's function is also explicitly made in Eq. (|54|) . allow- 
ing one to calculate, in principle, any physical observable. 

In Sec. lIVl we show that the problem simplifies greatly 
in the limit of negligible Landau level mixing. First, for 
locally flat potentials (away from saddle points or bottom 
of potential wells), we find that the mth Landau level 
acquires a dependence on the position R, according to 
the simple formula, 

e m ,±(R) = 5+(R) ± yWv^)2+ [fl-(R)] 2 . (5) 

Here v+(R) and u~(R) arc rcnormalizcd effective po- 
tentials that are simple functionals of the bare scalar 
and mass potentials ; for their definitions in terms of 
V s and V z , see Eq. (|64|) and the associated discussion in 
Sec. IIII Al Second, when curvature of the potential is in- 
cluded, wc find that simple analytic solutions for several 
parabolic models can be obtained. In particular, for cir- 
cular parabolic scalar potential V s (r) = (l/2)Uo(x 2 +y 2 ), 
the discrete energy spectrum (in terms of Landau level 
index m > 1 and an extra quantum number n, which is 
a positive integer > 1) reads: 

Vm + l 2 B U {m + n + 1/2) (6) 

(we have assumed fl c 3> l%Uo above). Apart from the 
well-known anomalous Landau-level quantization with 
respect to quantum number m, this result is very reminis- 
cent of Fock-Darwin states for standard 2DEGs with re- 
spect to the linear dependence in the integer n. More in- 
terestingly, for circular parabolic mass potential V z (r) = 
(l/2)Uo(x 2 + y 2 ), the discrete energy spectrum displays 
now an anharmonic form 

E m , n = ± yj (ffl^ + (l 2 B U [m + n + 1/2]) 2 

(7) 

that was not obtained to our knowledge. Generalization 
to non-circularly symmetric parabolic potentials is also 
readily obtained, as well as for the combination of uni- 
form scalar and parabolic mass terms (and vice versa), 
with detailed calculations appearing in several appen- 
dices. However, we show that potentials that combine 
sizeable spatial variations in both scalar and mass terms 
are in general not analytically tractable, even in the high 
magnetic field regime, except for the lowest Landau level. 



In Sec. El we make explicit the connection of our for- 
malism to the so-called deformation (or Weyl) quantiza- 
tion, which corresponds to the proper way of quantizing 
the dynamics in phase space. For two-dimensional prob- 
lems in a magnetic field, the vortex state formalism is in 
fact performing a mixed representation of phase space in 
terms of the two-dimensional coordinates of the center of 
mass together with a discrete quantum number associ- 
ated to Landau levels while standard Weyl quantization 
would introduce a four-dimensional description in terms 
of positions and momenta of the electron. This latter 
choice is however unpractical for the high magnetic field 
regime and this shows that the vortex states are most ro- 
bust in this regime. As should be expected, in the limit 
of infinite frequency (lu c — > oo for the 2DEG or Q c — > oo 
for graphene), Landau levels become fully decoupled, and 
the quantum dynamics reduces to a unidimcnsional one 
in terms of the two vortex coordinates, acting as conju- 
gate variables. An effective one-dimensional picture of 
motion is thus rigorously obtained, overcoming certain 
rcgularization problems of the path-integral technique^ 

Finally, in Sec. IVI1 we provide generic expressions for 
the LDoS in an arbitrary scalar or mass potential that 
can be described locally up to its first-order derivatives 
(generalized graphene drift states). Regarding recent ex- 
perimental findings, we show that: 

• positions, amplitudes, and widths of the LDoS 
peaks qualitatively depend on the dominant type 
(scalar or mass) of local potential, see, e.g., Eq. 

(d; 

• as the tip scans the surface, the LDoS peak en- 
ergy of the mth Landau level follows the effective 
potential given in Eq. (|5|) , see Fig. [TJ so that the 
resulting energy variations shrink with increasing 
m, in agreement with the experimental findings for 
graphene^ and standard 2DEGji 

• on the contrary, the width of the LDoS peaks at 
fixed tip position grows with increasing m (roughly 
as s/m), as observed in Rcf. |Tl| for graphene. 
Such a dependence is also expected for the ordi- 
nary 2DEG. 



II. FREE HAMILTONIAN : VORTEX STATES 
OF GRAPHENE 

A. Vortex states for the standard 2DEG 

Before investigating the case of graphene under mag- 
netic field, wc briefly recall the vortex states for the case 
of the non-relativistic 2DEG. This introduction will be 
useful to show that many physical and technical aspects 
of the 2DEG can be directly transposed to the case of 
graphene (studied in the next subsection). 

A single free electron of effective mass m* confined 
in a (xy) two-dimensional plane and subjected to a uni- 
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form magnetic field pointing in the perpendicular direc- 
tion B = Bz is described by the Hamiltonian 



H- 



2DEG 



n 2 

2m* 



K + ng 

2m* 



(8) 



Then, the eigenvalue problem H2deg^ = £^ leads to 
the well-known quantization of the kinetic energy into 
Landau levels, 



m + - ) 



(9) 



with the cyclotron pulsation lo c = \e\B/m*c — h/m*l 2 B 
and m > a positive integer (here / b = \/hc/\e\B is the 
magnetic length). It is important to note here the large 
degeneracy of the Landau energy levels e m . Indeed, for 
the motion of an electron in the two-dimensional plane, 
one expects at least two quantum numbers since there 
are two degrees of freedom. The degeneracy means that 
there is a great freedom in the choice of the second (de- 
generacy) quantum number, or equivalently, in the choice 
of a basis of eigenstates ^ . Consequently, there exist in 
the literature different ways to derive the energy quan- 
tization, Eq. ©. Eigenstates characterized by a pecu- 
liar symmetry of the (gauge-invariant) probability den- 
sity 1^1 2 are preferentially chosen in many contexts. For 
instance, the Landau states, with a conserved momentum 
as the degeneracy quantum number, are translationally 
invariant in one direction^ 2 - Circular eigenstates charac- 
terized by a rotation invariance around the origin^ are 
also well known and often used. It is worth stressing 
that the real difference between the Landau states and 
the circular states is not the gauge because both kinds 
of states can be obtained in any gauged The real dif- 
ference is in the choice of the gauge-invariant quantum 
numbers, which are intimately related to the symmetry 
of the probability density \^f\ 2 . 

Importantly, both Landau and circular eigenstates do 
not reflect the symmetry of the cyclotron motion around 
an arbitrary point R = {X, Y) in the (x, y) plane so that 
the consideration of the classical limit with these sets 
of states is rather tricky. Since they do not correspond 
to the classical picture of the motion, it is difficult to 
appreciate the wave-particle duality. By imposing that 
the probability density |>I'| 2 of the eigenstates has the 
same symmetry as the cyclotron motion, i.e., is a function 
of |r — R| only, we gei£Z- the so-called vortex states, given 
in the symmetrical gauge (A = Bz x r/2) by 



*m,R(r) 



x exp 



1 



x-X + i(y- Y) 

V%3 



lB\2-Km\ 
(x - Xf + (y - Y) 2 + 2i(yX - xY) 
472 



.(10) 



For practical convenience, we shall now use the Dirac 
bracket notation by writing ^ m ,n.{ v ) = (r|m, R). Eigen- 
states (fTU)) of Hamiltonian (|5J|, associated with energy 



quantization (|9]), are characterized by the set of quan- 
tum numbers |m, R), where m is a positive integer related 
to the quantization of the circulation around the vortex 
and R = {X, Y) is a continuous quantum number corre- 
sponding to the vortex location in the plane [note with 
Eq. (fTU| the "vortex" -like phase singularity at r = R 
for m > 1, which justifies the chosen denomination for 
the set of states] . These localized wave functions clearly 
encode the classical cyclotron motion around the guiding 
center R quantum mechanically. The vortex states form 
a semiorthogonal basis, with the overlap 



(mi,Ri|m 2 ,R 2 ) 



(RilR, 



(11) 



where 



(R1IR2 



exp 



(Ri - R 2 ) 2 - 2iz ■ (Ri x R 2 ) 



4/ 2 



(12) 

An important property is that the states (fTU)) present the 
coherent character with respect to the degeneracy quan- 
tum number R, i.e., they satisfy coherent states algebra. 
Note that these states are however eigenstates of the free 
Hamiltonian associated to the Landau-level index m, and 
form more precisely a scmicoherent basis with respect to 
the quantum numbers (m,R). In particular, they also 
obey the following completeness relation 



'rf 2 R 
2ttI% 



m=Q 



|m,R)(m,R| = 1. 



(13) 



According to this relation (| 13[) and general unicity prop- 
erties of the decomposition onto coherent states^ it 
is possible to expand arbitrary states or operators in 
the vortex state representation. Hence, despite being 
nonorthogonal, the set of states |m,R) with m > does 
form a basis of eigenstates, as the Landau and the circu- 
lar states. 

Besides providing a clear quantum mechanical dual of 
the classical cyclotron motion, there are several good 
reasons to prefer specifically the vortex states over an 
orthogonal set of eigenstates to study the process of lift- 
ing of the Landau level degeneracy in the presence of 
a smooth arbitrary potential. First, in contrast to the 
Landau states or circular states, the vortex states do not 
impose a symmetry to the degeneracy quantum number, 
and thus permit a great adaptability to the spatial vari- 
ations in the local electric fields, coming from either ran- 
dom impurity donors, confinement potentials, or macro- 
scopic voltage drops (in a nonequilibrium regime). This 
property leads to advantages in terms of computability 
since it is possible in the vortex representation to cal- 
culate and classify Landau-level mixing processes in a 
simple and natural manner (this will be illustrated in 
Sec. MI A[) . Second, at a more fundamental level, the 
vortex states are expected to be quite insensitive to any 
kind of smooth perturbations, since the quantum num- 
ber m has a purely topological origin in the vortex rep- 
resentation (for the Landau states or circular states, the 
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quantization of the kinetic energy comes either partially 
or entirely from the condition of vanishing of the wave 
function at infinity, what makes them much less robust 
to perturbations as a result of their nonlocality). Owing 
to this quantum robustness, the vortex states are thus 
naturally selected by the dynamics in the presence of a 
smooth potential with an arbitrary spatial dependence. 
They appear to be much more stable than their superpo- 
sitions (for instance, the Landau states) since they are the 
only states surviving under the action of such an interac- 
tion potential without any internal symmetry. Interest- 
ingly, the vortex states are also the best states to describe 
the transition from quantum to classical. Despite be- 
ing fully quantum, they thus encode de facto classicality 
properties and insensitivity to openness of the system. 
Therefore, they provide the best playground to under- 
stand the mechanisms of irreversibility, decoherence and 
dissipation in high magnetic fields. We will comment on 
this point in more detail later, in Sec. IIVDI 



B. Graphene vortex states 

We now come for good to graphene, which is described 
in the absence of potential by Hamiltonian (Q}. By 
searching the wave functions under the spinorial form 



* = 



with 



we get the following equations: 
(ti x - ifly^J w 



E 
— \ 

v F 

E 



(14) 

(15) 

(16) 
(17) 



with E the energy eigenvalue. Getting rid of the com- 
ponent u we get the Schr6dingcr-type equation for the 
component w, 



E 



(ll X + ilLyj (ll X - illy 

Using that 



— w. (18) 



n X) iL 



it) 



e\B 



'l 2 ' 

1 B 



we find that Eq. (|T8|) reads 

Il 2 w = Ew 

with 



E 



E 



I 2 ' 

1 B 



(19) 



(20) 



(21) 



By posing E ~ 2m* 'e in Eq. (|20)) . where e has the dimen- 
sion of an energy, we directly recognize the eigenproblem 
for a free 2DEG under magnetic fields discussed in the 
former section. This mapping shows that there is also a 
great freedom to choose a basis of eigenstates in the case 
of graphene. In the following, we introduce the analog of 
vortex states, Eq. (|10j) . for graphene. 
From Eq. @, we directly deduce that 



E = (2m + 1) h 2 /l\ 



(22) 



Therefore, we get that the energy eigenvalues of the 
graphene Hamiltonian are 



En 



A 



hV2 



Vf 



\\fmKl r 



(23) 



where A is a band index, which is equal to ±1 if m > 1, 
and if m = 0. We see that the energy levels are no more 
equidistant in energy and that the characteristic energy 
for graphene reads Hfl c cx \f~B~ instead of hoj c ex B for 
2DEGs. The component u of the spinorial wave function 
$ is straightforwardly obtained from the knowledge of 
the component w by using Eq. (|16[) . The corresponding 
normalized graphene vortex states are thus 



* m ,R,A(r) = 



1 



A* 



m-l,R 



(r) 



(24) 



Within the Dirac notation, the set of vortex quantum 
numbers we shall consider for graphene takes therefore 
the form 



|m, R, A) 



1_ _ / A|m-1,R) 



(25) 



The label A which characterizes the spinorial structure of 
the eigenvectors appears here as an additional quantum 
number with respect to the 2DEG. 

Using the semiorthogonality property, Eq. ([TT1) . of 
the vortex states, we can easily check that the graphene 
vortex states present the same property as their "non- 
relativistic" counterparts. Indeed, we have 



(toi,Ri, Ai|m 2 ,R 2 , A 2 



1 



1 



Ai(mi-l,Ri|\ _ / A a |m3-l,Ra) 
-i(mi,Ri| / 1 i|m 2 ,R 2 ) 



= 6 r , 



(Ri|R 2 



AiA 2 + 1 



yT+\xT\^T+\X2\ J 

= ^mi,m 2 (Rl|R2)<JAi,A 2 - (26) 

For convenience, in the next section we shall condense 
the full set of quantum numbers |m, R, A) into the single 
notation \i>). Therefore the sum over quantum numbers 
v will stand for 



E 



m=0 A 



(27) 
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It is finally straightforward to prove that the set of 
graphene vortex states \m, R, A) obeys a completeness 
relation, which reads 



EMM -/^ EE it 

1/3 m=0 A 



A 2 |m- l,R)(m - 1,R| -iX\m - l,R)(m,R| 



i\\m, R) (m — 1, R| 



|m, R) (m, R| 



d 2 R +°° 



E 



27T/2 

^ m=0 



1 
1 



|to,R)(to,R| 

|to,R)(to,R| 



(28) 



where we have used the completeness relation ([13]) satis- 
fied by the vortex states. 



III. GENERAL FORMALISM FOR A SMOOTH 
POTENTIAL 

A. Matrix elements of the potential 

In order to investigate the effect of a smooth potential 
under magnetic field, we shall naturally project the dif- 
ferent contributions of Hamiltonian in the graphene 
vortex representation. Although being basic, this projec- 
tion sheds already interesting light on the different pro- 
cesses at play and shows the essential differences between 
the different kinds of potentials that may be encountered 
in graphene, see Eq. (j4]). Using Eq. (|25l) . the matrix 
elements of the diagonal part of the potential (i.e., as- 
sociated to scalar and mass potentials) can be written 
as 



(^l|K 1 ia g k 2 )=[(l + |A 1 )|(l + |A 2 |)]- 1/2 

x {AiA 2 (mi - l,Ri]K + 14lm 2 
+(mi,Ri|V r ,-V r ,|m ai Ra)}. 



l,R 2 



(29) 



The off-diagonal terms of the potential (i.e., the random 
vector potential contribution) give rise to the following 
matrix elements: 



HKffh) =i[(l+ lAxIXl + IAa)]- 172 

x {Ai(toi - ljRilVa, - iV y |m 2 ,R 2 ) 
-A 2 (mi,Ri|T4 +iV v \m 2 - 1,R 2 )}. 



(30) 



We have shown in Ref. |27J that it is possible to evaluate 
exactly the matrix elements of a smooth function V(r) 
in the vortex representation [provided that V(x, y) is an 
analytic function of both x and y] and write them as a 
series in powers of the magnetic length I b , 



(m 1 ,R 1 |y|m 2 ,R 2 ) 



(R1IR2 



(R12) (31) 



with R i2 = [Ri + R 2 + i(R 2 - Ri) x z] /2 and 



,(R) 



d 2 V** mu n(v)*m 2 Ari)V(v) (32) 



(33) 



(n) 



y-^ ( TO 1 + ^mi+fc,m 2 +3-fc 



k=0 



Vmi!m 2 ! k\(j-k)\ 



x (d x + id Y f (dx - id Y ) 3 ~ k V(R 



Clearly, the use of an analytical expansion around the 
complex point Ri 2 in Eq. (|31[) puts some constraints 
on the types of potential that can be considered in the 
present formalism. We emphasize that relation (|3"Tj) 
holds for any physical potentials V (which are neces- 
sarily smooth functions of the space variables). In 
contrast, pointlike (i.e., zero-range) potentials involving 
Dirac delta functions which represent toy models simu- 
lating short-range potentials can not be treated within 
the present formalism. If the magnetic length Ib corre- 
sponds to the shortest length scale [here, basically, Ib 
has to be compared with the characteristic length scale 
of spatial variations in the function V(R), see Eqs. (|33|) 
and (|34p]. we see that we have naturally ordered the dif- 
ferent contributions to the matrix elements by their order 
of magnitude in high magnetic fields. 

At leading order (l B -> 0), we get from Eqs. ([29 ]l -([3"i ]) 
for coinciding vortex positions Ri = R 2 = R, 

(^llKiiagl^) « 5 mu m 2 [S\ u x 3 Vs (R) - Sx!,-X 2 Vz (R)] • 

(35) 

We remark that in the limit Ib — > the diagonal ele- 
ments V s and V z of V do not introduce a mixing between 
Landau levels. For smooth functions V x and V y we get 
in the same limit Ib — > 0, 



HKrfih) «- i [(l + |A 1 |)(l + |A 2 |)]- 1/2 
x {X 2 S mum2 -i [V X (R) + iV y (R)] 

-Ai5 mi _i, m2 [V X (R) - iV y (R)}} . (36) 

We note with Eq. (f36|) that the off-diagonal elements V x 
and Vy do mix adjacent Landau levels already at lead- 
ing order in Ib, in contrast to the diagonal elements V s 
and V z of V. This difference clearly calls for a differ- 
ent treatment of the diagonal and off-diagonal parts of 
the total potential V. Off-diagonal contributions can be 
treated perturbatively at high magnetic field by assum- 
ing that V x and V y are small in amplitude in addition of 
being smooth functions at the scale Ib- Such a constraint 
on the amplitude can be relaxed in the treatment of the 
diagonal contributions of V. 

The next (sub-dominant) contributions of order Ib to 
the matrix elements of Vdiag are proportional to 

(d x + idy) [V S (R) - V Z (R)} 

+AiA 2 ^r[^(R) + V Z {R)]} + c.c.(l O 2), (37) 
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where the notation c.c.(l <H> 2) means taking the complex 
conjugate and exchanging the indexes 1 and 2 of the for- 
mer expression. This contribution induces a mixing be- 
tween both adjacent Landau levels and band indices A. 
Moreover, the mixture of positive- and negative-energy 
components stems from both components V s and V z of 
the potential energy. It is interesting to note that for 
a large Landau-level index, the mixture arising purely 
from V s (i.e., taking the mass term V z = 0) gets neg- 
ligible when AiA 2 = — 1. For instance, when mj and 
m 2 3> 1, we have 



%/toi+T + AiA 2 i 



fm{ (1 + AiA 2 ) 



(38) 



for the component <5 mi +i. m2 of the matrix elements, Eq. 
(|3"T|) . associated with V s . On the other hand, the band 
mixing becomes significant for mi and to 2 close to 0. 
Specific signatures resulting from this interband mixing, 
such as Zitterbewegung (or trembling motion) in a mag- 
netic field, have been discussed in the literature^ - — By 
looking at next-order contributions in l\ for the matrix 
elements, we note that interband mixing occurs also with 
the second derivatives of a pure scalar potential V s with- 
out mixing the Landau levels. These mixing processes 
will be analyzed further in Sec. IIVI 



B. Green's- function formalism 

The nonorthogonality of the graphene vortex states 
preventing us to build a wave-function perturbation the- 
ory, we shall instead use a Green's-function formalism to 
get a more quantitative insi ght on the effect of a smooth 
potential, following Refs. |28| and [26|. Although the 
derivation of the equations of motion for the graphene 
Green's function is very similar to that for the 2DEG 
Green's function, we shall nevertheless describe the prin- 
cipal steps with some detail here, in order to make this 
paper self-contained (we shall, however, not reproduce 
the very technical details). 

Retarded and advanced Green's functions are, respec- 
tively, defined as 

G R (x 1 ,x 2 ) = -i0(t 1 -t 2 )({ij(x 1 ),^(x 2 )}), (39) 

G A (x 1 ,x 2 ) = i6(t 2 -h){{il>{xr),il>\x 2 )}), (40) 

where {, } means the anti-commutator, and the Hcavi- 
side step function [i.e., 6(t) = for t < and 6(t) = 1 for 
t > 0] . The averages are evaluated in the grand canonical 
ensemble. The Green's functions relate the field opera- 
tor ip( x i) °f the particle at one point x\ = (ri,ti) in 
space-time to the conjugate field operator 7/^(2^) at an- 
other point x 2 = (r 2 ,t 2 ). The field operators tp(x\) and 
ipi(x 2 ) are expressed in terms of the eigenfunctions ^(r) 
and eigenvalues E v as 



,iE v t 2 /h 



(41) 
(42) 



where c\, and c v are, respectively, the creation and de- 
struction operators. 

As a basis of states, we shall then use the graphene 
vortex states \v) = \m, R, A) which are eigenstates of 
Hamiltonian Ho [Eq. (JTJ]. It is worth noting that, 
although these states \v) are nonorthogonal, the as- 
sociated creation and destruction operators cj, and c v 
obey the usual algebra with the anti-commutation rules 

{4i)4a} = {^1^2} = and {^i>4 2 } = ^in- 
completeness relation (f2"5)) allows us to express 

the Green's function in the graphene vortex rep- 
resentation, which we note G R ' A (v\, 1\) v 2 , t 2 ) = 
G H,j4 (mi,Ri, Ai,ti;TO2, R 2 , A 2 , i 2 ). Transposing its def- 
inition originally made in terms of the electronic coordi- 
nates (r,t) into the vortex language, the latter Green's 
function gives the probability amplitude for a vortex with 
circulation mi and band index Ai that is initially at posi- 
tion Ri at time t\ to be at point R 2 at time t 2 with a new 
circulation m 2 and a band index A 2 . After Fourier trans- 
formation with respect to the time difference t = t\ —t 2 , 
the Green's function (denoted by Go) corresponding to 
Hamiltonian Hq [i.e., Hamiltonian with V = 0] are 
written in the energy (cj) representation as 



rtR.A, \ ^mi,m 2 <^Ai,A 2 (Rl|R2 

G [y Y \v 2 ) = - 



E 



mi ,Ai 



(43) 



Retarded and advanced Green's function in the pres- 
ence of the smooth potential V are obtained from Dyson 
equation, which takes the following form in the v repre- 
sentation (we again considered the Fourier transform of 
Green's function with respect to time difference) 

(w - E miM ± i0 + ) G R ' A { Vl ; va) = (vi\v 2 ) 



Y,V„ u „ 3 G R ' A (v s ;v 2 ). (44) 



Here the general matrix elements V UlW2 = {vi\V\v 2 ) = 
(Ri|R 2 ) fm 1 ,A 1 ;m2,A 2 (Ri2) are given by expressions 
(p4| . For 7^0, the graphene vortex Green's function is 
generally no more diagonal with respect to the quantum 
numbers m and A, and the mixing between the different 
quantum numbers depend on the characteristic proper- 
ties of the potential V. However, it turns out that, as 
a result of the coherent states character with respect to 
vortex position R encompassed within overlap (|12[) . the 
propagation of the graphene vortex Green's function with 
respect to vortex positions Ri and R 2 is constrained to 
necessarily take the form 



G{V\] V 2 ) — (Ri|R 2 ) <7rai,A i; ra2,A 2 (R12) , 



(45) 



similarly to the matrix elements of the potential [see Eq. 
(|31j)]. Such exact dependence, Eq. (|45|) . can be de- 
rived from Dyson Eq. (|44[) in the same way as done 
in Ref. H|| Remarkably, it implies that the nonlocal 
graphene Green's function G{v\]V 2 ) will be entirely de- 
termined once it is known at coinciding vortex positions 
Ri = R 2 = R, and this result holds irrespective of the 
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potential V. It is then sufficient to consider Eq. (|44j) 
for coinciding vortex positions. Because the derivation is 
the same as for the 2DEG, we briefly outline here the last 
step leading to the final equation of motion governing the 
function Om li A 1 ;m 2 ,A2 an d refer the reader to Sec. II 
of Ref. |23 for the mathematical details. The nonlocal 
dependencies of the functions G(ys; v-i) and V Vl - U3 on the 
vortex positions which are known according to relations 
(|45j) and (|31l) arc exploited to evaluate the integral over 
the continuous variable R3 on the right-hand side of Eq. 
([44]) . This integral then transforms into a series expan- 
sion in powers of Is- We obtain that Dyson equation 
for the retarded graphene vortex Green's function g(R) 
(from now on, we drop the R upperscript associated to 
retarded) corresponding to Hamiltonian ([3]) reads 

(U) — £ mi| \, + «0 + ).g mi ,Ai;m 2 ,A 2 (R) = <5mi,m 2 ^Ai.A 2 

+ °° / I \ 2k 1 

fc=0 w^/ K - m3 . M 

x (d x +«9y) fc g m3 ,A 2 ;m 2 ,A 2 (R) (46) 

with E m x = ^E m = \h\Jlmvp Hq. 

Another important aspect of the change in function 
, which appears clearly with the form (|46p of Dyson 
equation and with expressions (|33[) and (|34[) for the ma- 
trix elements of the potential taken at coinciding vor- 
tex positions, is that the nonanalytic dependence of the 
nonlocal graphene vortex Green's function G{y\\v-i) on 
the magnetic Ib has been entirely extracted [in formula 
(|45p. this nonanalytic dependence is only contained in 
the overlap (R1IR2)]. In other terms, the function g(R) 
is obviously analytic in Ib and thus well behaves in the 
scmiclassical limit of zero magnetic length (Ib —> 0). This 
property can be used to solve Eq. (|46[) order by order in 
powers of Ib and thus to provide a semiclassical expan- 
sion of the graphene vortex Green's function g r n, 1 ,\ 1 ;m 2 ,\2 
as 



?mi,Ai;m 2 ,A 2 



+00 

E 

3=0 



l_B_ 

s/2 



"mi ,Ai ;rri2 ,A2 ' 



(47) 



Because the series, Eq. ((47)) . is then only asymptotic in 
nature (the obtained solution holds in the limit Is — > 0, 
but is not controlled at finite Ib), we aim here at solving 
directly and non-perturbatively in Is Dyson Eq. (|4l)l) . 

For this purpose, we have found in Ref. Hfj] that it is 
very convenient to introduce the simultaneous changes in 
functions, 

.9 mi ,A i; ™ 2 ,A 2 (R) = e-^^ AR g mitXl . im2t x 2 (R), (48) 
« mi ,A i;m2 ,A 2 (R) = e- (i « /4)AR v mi , Ai;m2 ,A 2 (R), (49) 



where the symbol Ar means the Laplacian operator 
taken with respect to the vortex position R. After 
substitution of these expressions (|48[) and (|49| into Eq. 
(|46p . we get a new equation for the unknown function 
.9m i , Ai ;m 2 ,a 2 (R) with a higher-order differential operator 
than the one appearing on the right-hand side of Eq. (|46|) 



(tJ — E mii \ t + «0 + )g mi ,Ai;m 2 ,A 2 (R) — ^Ai,A a ^mi, m 2 

5m3,A 3 ;m 2 ,A 2 

(R), (50) 



where the symbol * stands for the bi-diffcrcntial operator 
defined by 



exp 



4 (** 



>Y - ^Y~3 X 



(51) 



The arrow above the partial derivatives indicates to 
which side the derivative acts. Note that the passage 
from Eq. ((46]) to Eq. (|50|) is more straightforward by 
going to Fourier space (see Appendix A of Ref. 1261 ) . It is 
worth mentioning that, by starting from the other Dyson 
equation (i.e., formally G = Gq + GVGq) and following 
the same steps as detailed previously, we can derive a 
second equation satisfied by the function g, 

((J — E m2t x 2 + 20 + )(7 mii Ai;m 2 ,A 2 (R) = Sx 1 .X 2 Sm 1 ,m 2 

(R). (52) 

«13,A 3 

The particular forrn^ of exact Eqs. ([50|) - ([52)l . reminis- 
cent of the so-called star-product, will be further used 
and commented in Sees. IIVI and fVl 



In order to compute local physical obscrvablcs such as 
the local density of states, we need to express Green's 
function in terms of the electronic positions r. The elec- 
tronic Green's function is a 2 x 2 matrix in the pseudo- 
spin space, and is defined as G(r,r') = (r|G|r'). At a 
practical level, it is useful to directly relate the nonlocal 
electronic Green's function to the local graphene vortex 
Green's function <?mi,Ai;m 2 ,A 2 (R) (at coinciding vortex 
positions) or alternatively to the modified vortex Green's 
function <? mij Ai;m 2 ,A 2 (R)- First, the electronic Green's 
function can be straightforwardly linked to the nonlo- 
cal graphene vortex Green's function G(y\\v-i) through 
a change in representation which is performed by using 
twice completeness relation (|28|). Then, using Eq. (|4"5)) 
and following the calculations made in Ref. [28| for the 
2DEG, we get the following relation 
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G(r,r» 



rf 2 R 

nl% 

D mi,Al 1712-^2 



E E 



AiA 2 ^ 2 _ ljR (r')* mi -i, R (r) -iAi^ 2iR (r')* mi -i,R(r) 
i^* _i, R (r')* mi ,R(r) ^ J , R (r')*m I ,R(r) 



xe 



^mi , Ai ;m2 ,A2 (^V) 



(53) 



where the functions \E , m ,R(r) correspond to the so-called vortex wave functions written in Eq. (JTUJ). Inverting 

expression (|48p . i.e., writing g mi! > i;m2 ,A 2 (R) = e^ B / 4 ^ AR g mi .x 1 ;m.2,x 2 ( R ) an d inserting this result into Eq. (|53p . we 
get after integrations by parts (so that the operator involving the Laplacian acts on the product of wave functions 
rather on the local vortex Green's function) 



G(r,r» = 



d 2 K 



2nl 



rE E 



A R ( AiA 2 ** _ liR (r')* mi -i )R (r) !R (r')* mi -i, R (r) 



mi,Ai m 2 ,A 2 



»A 2 *,* i ,_ 1 , R (r')*m 1 3(r) 



*™ 2 ,R(r')* mi ,R(r) 

9mi,Ai;m2,A2 ( R ) 



(54) 



Because the functions gm. 1 ,x 1 :m 2 ,X2 ma y depend on Ai and A 2 , the electronic Green's function G(r, r') possesses, in 

general, off-diagonal elements. The above equation is a central one, because it shows that any physical observable 

can be computed from the knowledge of the local vortex Gr een's function ffym.Ai;m 2 .A2 (R). 

I — ' 



IV. HIGH MAGNETIC FIELD REGIME 
A. Regime of negligible Landau- level mixing 

While Eqs. (|50| - (|52)l can, in principle, be consid- 
ered for any magnetic fields, we shall investigate here 
the regime of high magnetic field only, for which Landau 
level mixing can be safely neglected. This regime can be 
reached under reasonable conditions (i.e., for fields on the 
order of 1 T or higher) provided that the potential land- 
scape is sufficiently smooth. Indeed, Landau-level mixing 
processes are described within Eq. ([50)1 by the matrix el- 
ements Vmi,Xi;m 3 ,X 3 with mi ^ 7713. From the expressions 
of the matrix elements of the potential coupling adjacent 
Landau levels calculated in the vortex representation in 
Sec. IIII AI we can formulate a clear quantitative criterion 
for neglecting Landau-level mixing due to the diagonal 
contributions of the potential V in graphene, 

/ B |V R V(R)| < (Vm + 1- Vm)hn c . (55) 

In graphene and for a field of 5 T, we have Kl c = 
H\/2vf /Ib ~ 80 meV and Ib — 11 nm. Recent exper- 
imental STS measurements of the spatial dispersion of 
Landau levels in epitaxial graphene^ give at most typ- 
ical linear variations in 5V ~ 5 meV on length scales 
£ ~ 20 nm. Thus l B \W\/Ml c < l B 5V/^hfl c ~ 0.03, a 
very small number indeed, so that the limit of negligible 
Landau level mixing is well obeyed. We shall furthermore 
suppose that the Landau level mixing processes due to 
the off-diagonal part of V are small. According to Eq. 
(ptrJ)) . this implies 

\V Xt y(B.)\ < (v^TT - y/m) tin c . (56) 



Under inequalities (|55p and (|56[) . Landau- level mixing 
processes due to the spatial variations in the scalar po- 
tential V s and of the random mass V z or to the spatial 
fluctuations {V x and V y ) of the vector potential are small 
and can be accounted for perturbatively on the basis of 
Eq. 0. 

Henceforth, we shall concentrate on the main relevant 
processes occurring at high magnetic field in a smooth 
potential. In this regime, the Landau-level degeneracy is 
principally lifted by the presence of both the potentials 
V s and Vz, which give rise for to > 1 to the following 
diagonal (mi = to 2 = to) matrix elements in the vortex 
representation 

^m; Ai ; A2 (Jty — ^mi ,m 2 ^mi , Ai \m, 2 , A2 ( R ) 

= <*Ai,A 2 «m( R ) + <*Ai -A 2 V~(R), ( 57 ) 

where the diagonal and off-diagonal components of the 
potential matrix elements in pscudospin space, respec- 
tively, read 

vt(R) = \J d 2 V [\V m Mv)\ 2 (V s (v)-V z (v)) 

±\* m - hn (r,)\ 2 (V s (r,) + V z (r,))] (58) 

-5g^(S) i {"-w-"-w 

±J?-(V S (R) + V Z (R))}. (59) 

To write down expressions (J57j) - (J59]) , we have used Eqs. 
(|2T)1) and (|3"2"j) - (|3"^|) . We notice that even a scalar poten- 
tial V s introduces a coupling between the bands A = ± 
for to > 1 through its nonlocal differential contributions 
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arising with j > 1. For instance, a quadratic scalar po- 
tential generically mixes the positive and negative energy 
components, even in the absence of a mass term (V z = 0). 
The case to = has to be treated as a special case since 
there is only one band (A = necessarily). The matrix 
elements for the lowest Landau level m = read 



d 2 V ^o,u(v)\ 2 [Vs(r])~V z ( V )} (60) 



= E^ff Ar)V,(R)-^(R)]. (61) 

We have seen previously that Dyson Eq. ([50]) is greatly 
simplified when considering modified matrix elements 



^m, Ai ;m, A2 (R) 



-0i/4)Ar 



Wm ) Ai;m 1 Aa(R)) which consti- 



tute the effective •potential in Landau level m. Using re- 
sults given in the Appendix B of Ref. H(| we get the 
action of the exponential differential operator onto the 
product of two vortex functions with identical Landau 
level m and positions r: 



K m (n 



- (i - /4)AR |* m ,R(r) 



-A s (R,- 



2 nl 



irm\l 2 B ds m 



1 + s 



(62) 
,(63) 



s=0 



with A s = (1 — s)/(l + s). Thus, the diagonal and off- 
diagonal effective potentials (in pscudospin space) read 
for to > 1, 

S±(R) = ^[^(R-ffWfaJ-V,^)) 

±lf TO _i(R-i|)(V;(t/)+V,(^))]. (64) 



We emphasize that formula (|64j) is non-perturbative in Z b 
and possibly applies for potentials V s and V z with sizeable 
variations at the scale of Ib- The effective potential in 
the lowest Landau level is also readily obtained as: 

60 (R) - J d 2 V K (H - 77) [VM - V Z ( V )] ■ (65) 

Obviously, we find that the modified Green's function 
becomes also diagonal with respect to the Landau-level 
quantum number at large magnetic field (yet at finite 
magnetic field), 

gm 1 ,A 1 ;m 2 ,A 2 (R) = 5mi,TOa Smi,Ai,A 2 (R)j (66) 

and is determined for m > 1 by Dyson equation, 
(u - E mM + iO + ).g m; A i; A 2 (R) = <^Ai ,A 2 

+«m(R)*ffm;Ai;A 2 (R) + "m( R ) * ffm;-Ai;A a ( R ) (67) 



and for to = by 



(c^0+).g o (R) = l + u (R)*5o(R) 



(68) 



The other Dyson Eq. (|52|) generates the different equa- 
tion, 

(u - E mM + iO + ).g m; A i; A 2 (R) = <^Ai ,A 2 

+fl m ;A 1 ;A 3 (R)*U+(R)+5 m ;A 1 ;_A 3 (R)*«~(R) (69) 

for to > 1, and 

{oj + iO+)g (R) = l+g (R)*S (R) (70) 
for the lowest Landau level m = 0. 

B. Locally flat potentials 

Now, we aim at solving Eqs. (|6"T|) - (fT0"|) at leading or- 
der, which is vindicated when the potential is locally flat, 
i.e., when potential curvature is small. This calculation 
includes the case of one-dimensional potentials (i.e., glob- 
ally flat potentials), for which the solution presented be- 
low is exact. Indeed, as is clear from its explicit expres- 
sion (|51|) . the ★-bidiffercntial operator involves deriva- 
tives in two orthogonal positions. In case where the po- 
tentials U S (R) and V^(R) are purely one-dimensional po- 
tentials depending on the same coordinate, the function 
<?(R) will also only depend on the same and unique vari- 
able, so that the * product between the functions v and 
g reduces to the standard product of functions. In case 
of arbitrary spatial varying two-dimensional potentials, 
this constitutes a good approximation as long as temper- 
ature is higher than the energy scales associated to local 
curvature terms, see Sec. IIVDI for a general discussion. 
Dyson equation then is trivially solved, as the system of 
differential Eqs. (|67| -([70 |) transforms into a system of 
purely algebraic equations. Taking the difference of Eqs. 
(p3?|) and (|69p . we get for to > 1 the relations between 
the different components of g, 



<7m;+;- (R) = 5m;-; + (R) 



(71) 



g m; -,-(R) = 3 m;+;+ (R)-2^g-5 m; _ i+ (R). (72) 

v m (R) 

After simple algebra, we directly obtain the solution 

5m;Ai;A2 (R) = [w _ ^ +(R) + . Q+] [u _ em> _ (R) + i0+] 

x { [w - v+(R) + E mAl ] 6 XlM + v-(R)5 Xl ,-\ 3 } (73) 

with the poles (corresponding to the renormalized Lan- 
dau levels) giving the effective energies, 



£ m ,±(R) = S+(R) ± y/Efr + [i m (R)]. (74) 

For to = 0, the Green's function is characterized by a 
single pole and reads 



5o (R) 



1 



w-£o(R) + *o+' 



(75) 
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where £o(R) = ^o(R). 

Equation (fM|) . with the explicit expression for the 
renormalized potentials given in Eqs. (|64|) . provides the 
leading result for the local Landau-level energy in arbi- 
trary potentials of diagonal type (i.e., scalar or mass- 
like). This expression of course includes the case of a 
purely unidimcnsional (i.e., globally flat) potential as an 
exact particular solution, but is a very good approxima- 
tion for smooth disordered potentials, which can be used 
to analyze experimental STS results, as we discuss in 
Sec. EH 



C. Locally curved potentials 

This section presents the resolution of Dyson equation 
at next to leading order, by extending the above calcu- 
lation of the local vortex Green's function to the incor- 
poration of the effects of geometrical curvature in the 
potential landscape. It has therefore a two-fold purpose. 
First, it provides a crucial refinement of the previous ex- 
pression (|73|). that includes important quantum effects 
such as quantization of energy levels or tunneling as- 
sociated to the potentials V s and V z , which are clearly 
missed in the leading order guiding center Green's func- 
tion. Smaller energy scales associated to these physi- 
cal processes are now accessible, and the final expression 
will apply to arbitrary smooth potentials that are locally 
curved. These important aspects are discussed in more 
detail in Sec. IIVDI Second, in the special case of purely 
quadratic potentials (which thus have a global constant 
curvature), the calculation provides essentially the ex- 
act Green's function, from which one can gain interest- 
ing insights on the physics of confinement or tunneling 
in graphene. We thus obtain analytically the quantiza- 
tion spectra of parabolic quantum dots and show that 
the structure of energy levels qualitatively depends on 
the type of confinement (electrostatic or mass type). We 
henceforth assume that the diagonal potentials are locally 
well described up to their second-order spatial deriva- 
tives. 



1. Lowest Landau level: Solution with both curved scalar 
and mass potentials 

We start by considering the lowest Landau level m = 0, 
which is the simplest case to solve, as band indices are not 
involved. In that situation both locally curved V^(R) and 
V~(R) can be solved altogether (this is not the case for 
higher m > 1 states, as will be discussed in the next para- 
graphs). Actually, Dyson Eq. for the lowest Landau 
level is formally equivalent to the equations obtained^ 
for the 2DEG, as the electrostatic potential V(R) for 
the 2DEG is just formally replaced by the combination 
T4(R) ^ Vz(R) for graphene. Working in the next to 
leading order, i.e., keeping local curvature terms of order 
lg in the *-bidiffercntial operator, Eq. (|5lT) . we can di- 



rectly transpose the solution of Ref. [26| to the graphene 
case (for the method, see also Appendix [B] of the present 
paper), which reads 



5o (R) 



with 



+ oo i[^ (R)/ 7o (R)][t-T (t)] 



i I dt- 

10 



ro(t) 



;t[w-£o(R)+io 4 



cos(y / 7o(R)i) 



v/toTR 



tan(V7o(R)*)- 



(76) 



(77) 



The parameters 70 (R) and 770 (R) in Eqs. (|76|) and ([77)) 
are geometric coefficients characterizing the local effec- 
tive potential landscape «o(R) in the lowest Landau level: 



7o (R) 

%(R) 



R ' 



(78) 



'-f [{d 2 x V ){d Y V ) - (d X d Y V ) 2 

l -f [(d 2 x v, )(d Y v ) 2 + (d Y v a )(d x v f 

-2(d x dYV )(d x v )(dYV )] R . (79) 



The coefficient 70 (R) is directly proportional to the 
Gaussian curvature of the surface defined in the three- 
dimensional "space" XYZ by the equation Z = 
Vq(X,Y). Its sign reflects the local topology of the ef- 
fective potential: 70 (R) > indicates a locally elliptic 
potential with the presence of a local extremum (max- 
imum or minimum), while 70 (R) < corresponds to a 
locally hyperbolic (or saddle-shaped) potential. At the 
borders between the regions with curvatures with oppo- 
site signs, the potential is locally parabolic (the lines 
where the Gaussian curvature is zero are consequently 
called parabolic lines). For a complex disordered effec- 
tive potential landscape, one expects that surface regions 
with positive and negative Gaussian curvature alternate. 
Note that both cosine and tangent trigonometric func- 
tions in Eqs. (|76[) and (|77[) transform into their hyper- 
bolic counterparts in the case 70 (R) < 0. Equation ([To) 
thus provides a general approximation scheme in the low- 
est Landau level in the presence of arbitrary scalar and 
mass potentials that are locally well described by local 
curvature coefficients, Eqs. ([75)1 and (f79"j) . 

Now, in the particular case of purely quadratic scalar 
and mass potentials, i.e., V^(R) — V Z (R) = V s o — V z q + 
i [(R — Ro) • VrJ 2 (V s — V z ), with R chosen as the sin- 
gle point where the potential gradient vanishes, expres- 
sion (|76p yields the exact Green's function of the problem. 
In that situation, the parameter 70 in Eq. (|78[) becomes 
R-independent , 



To 



■f [d 2 x (V s - V z )d 2 Y {V s - V z ) - [dxMVs 



K)] 2 }, 

(80) 

and describes the uniform (global) curvature of the po- 
tential, while the R-indcpcndcnt part of the effective po- 
tential results from the simple relation u = wo(Ro) = 
fio(R) - 7?o(R)/7o = V s0 - V z0 + (i|/2)A R (V. - V z ). For 
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a confining potential, i.e., when 70 > 0, To(t) is a 2tt/^/q 
periodic function of time t. Direct Fourier analysis of ex- 
pression (|76p using the above relations can be done and 
shows (see Appendix [XJ that the entire energy spectrum 
necessarily decomposes onto discrete modes: 



151 The solution reads 



E , n = % + sgn(?7o)\/7o(2n + 1) 



(81) 



with n > a positive integer, yielding a harmonic- 
oscillatorlike spectrum for the parabolic quantum-dot 
model (in the large magnetic field regime considered 
here). The general form of this spectrum will be dis- 
cussed in the next section. In contrast, for 70 < 0, 
the vortex Green's function expressed in the time rep- 
resentation is no more periodic but decays on a time 
scale 1 / ^/— 70, due to the cutoff function 1 / cosh( % /— Tot) . 
These lifetime effects associated to negative Gaussian 
curvature are clear manifestations of quantum tunnel- 
ing in saddle-point potentials, and will be considered in 
a future publication where transport properties in high 
magnetic field will be considered. 



2. Arbitrary Landau level: Solution for a curved scalar 
potential combined with a flat mass potential 



For m > 1, the structure of Dyson Eq. (|67j) for 
graphenc differs from that for the 2DEG case because 
of the possible coupling between positive- and negative- 
energy bands. Two kinds of processes are actually at 
work here. First, non-zero mass potential V z directly 
couples the two bands, as is clearly seen from the lead- 
ing order Green's function in Eq. (|73l) . Second and 
less obviously, higher-order scalar processes can also in- 
duce band-mixing. Indeed, the effective off-diagonal po- 
tential in Eq. ([55)) reads in the small ^-expansion: 

v~(R) = -V z (R)-mf A R F z (R)+f A R F s (R)+0(Z|). 
Thus, even for an identically zero mass term (V z = 0), 
positive-and negative-energy bands are necessarily cou- 
pled by the second derivatives of the scalar potential. 

For reasons mentioned previously, one cannot analyti- 
cally progress for the m > 1 Landau levels in case where 
both scalar and mass potentials are strongly spatially 
dependent. In this section we therefore assume that 
the scalar potential varies in space with sizeable local 
parabolic dispersion, while the mass potential has much 
smoother spatial variations, so that local derivatives of 
the mass term are associated to tiny energy scales (the 
reversed situation, where the mass potential variations 
dominate the ones of the scalar potential, is considered 
below in Sec. lIV C~3|) . Since the calculation leading to the 
Green's function for graphenc is largely inspired from the 
2DEG's derivation^ details are produced in Appendix 



ffm;A i; A 2 (R) = ~ j di- 



ces (y 7 m(R)t) 

^ git(w-fm,e(R)-MO+) 



x[(l + eA 1 a m (R))5 AliA2 + e ( S m (R)5_ Al> A 2 ] (82) 

where the effective energy £ mi ±(R) is given by Eq. (|?4"|) . 
and 



a m (R) = 

MR) = 



E,, 



El + [t£(R) 

MR) 

El + [t£(R) 



(83) 
(84) 



The geometric parameters 7+(R) and 77+ (R) have the 
same definitions as in Eqs. ([78)1 and ([75)l . where «o(R) 
is simply replaced by the effective potential u+(R). The 
function r+ (t) has also a similar definition as in Eq. (f77|) 
now in terms of 7+(R). Again, the above expression (|B2")| 
is quite general, and can be used to describe arbitrary 
disordered (yet smooth) scalar potentials. A mass con- 
tribution may be present, but only with negligible spatial 
variations for the approximation to be valid. 

Now, in the particular case where the bare scalar po- 
tential is globally quadratic (i.e., has uniform curva- 
ture) and the mass potential is globally uniform, this 
expression provides the exact Green's function. A pos- 
sible parametrization of such potentials reads V^(R) = 
+ 5 [(R — Ro) ' Vr] 2 V s (with R chosen as the 
point where the scalar potential gradient vanishes) and 
V Z (R) = V z o- The Gaussian curvature of the scalar po- 
tential becomes then constant and independent of m, 



7+ = f [#xVJ*V. 



(d X d Y V s f 



(85) 



while the R-independent parts of the effective poten- 
tials read u+ = u+(Ro) = V s o + m(l%/2)AnV s and 
u~ = — Vzo + (1%/4)ArV s . Fourier analysis as done in 
Appendix [S] provides a spectrum of purely discrete en- 
ergy levels in the presence of 2D-parabolic scalar poten- 
tial, 



E n 



^JeI + (^)2+sgn( ?? +) V7+(2n+l). (86) 



This form of quantization is quite reminiscent of the 
Fock-Darwin spectrum for the non-relativistic 2DEG: be- 
sides the renormalization of Landau levels (labeled by 
the integer m) due to the R-independent part of the 
effective potentials v+ and the linear dependence 
in the second discrete number n provides an additional 
harmonic-oscillatorlike contribution. As a specific illus- 
tration for the case of a circular parabolic scalar potential 
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V s (r) = V s o + (l/2)U~o(x 2 + y 2 ) together with a zero mass 
term, one gets the following energy spectrum: 

E m , n = V s0 +l%U (m+n+±)±^ (M2 cV ^) 2 + {l 2 B U /2)\ 

(87) 

that we have already quoted in Eq. (|6|) in the large Q c 
limit. 



3. Arbitrary Landau level: Solution for a flat scalar 
potential combined with a curved mass potential 

We now consider the alternative solvable case of locally 
fiat scalar potential, together with a spatially dependent 
mass potential that can be locally well described by a 
quadratic expansion. Solution of Dyson Eq. (|67[) can then 
similarly be achieved, leading to the Green's function for 
m > 1, 



„i[K m (R)+iO + ]s 

5m;Ai;A 2 (R) = ~* / ds , 

COS ( W7m(R)s 



I cos[0 m (s)] 



w-t5+(R) + Ai^ 



■Sxi,Xa + ism[0 m (s)]8x lt -xA 



K m (R) 

with 6 m (s) = [r?-(R)/7"(R)][r-(s) - s] - sv~(R) and 
K m (R) = sgn [w - S+(R)] (w - «+(R)) 2 - ^ V2 (89) 



if \u - v+(R)\ > E m , and 



K m (R)=i [ W -5+(R)] 2 -El 



1/2 



(90) 



if |w — w+(R)| < E m . Details for the derivation of result 
([88} can be found in Appendix [C] The geometric param- 
eters 7" (R) and ^(R) in formula (|88|) are defined as in 
Eqs. ((7HJ) and JTSJ with u (R) replaced by -v~ (R). Tak- 
ing the imaginary part of expression (|88|) . we get that the 
local density of states vanishes when \lu — «+(R)| < E m , 
meaning that there are no states within this energy in- 
terval. 

Now, in the particular case where the bare scalar po- 
tential is globally uniform and the mass potential is glob- 
ally quadratic, this expression provides the exact Green's 
function in the absence of Landau-level mixing. A possi- 
ble paramctrization of such potentials reads V s (R) = V^o 
and V Z (R) =V M + \ [(R - Rq) ■ V R f V z (with R cho- 
sen as the point where the mass potential gradient van- 
ishes). The Gaussian curvature of the mass potential 
becomes then constant, 

7" = l f [d 2 x V z dlV z - {d x d Y V z ) 2 ] (91) 

while the R-independent parts of the effective potentials 
read 6+ = V r s0 -(^|/4)A R V Z and v~ = v~(Ro) = -V z0 + 



m(l 2 3 /2)AnV z . Fourier analysis as done before implies 
that the eigenenergies u> = E m n are determined by the 
implicit equation K m — v~ n = sgn(77^ ') (2n + 1) [we 
remind that the dependence on u is contained in K m , 
see Eq. [53] , leading to the following discrete energy level 
spectrum in the presence of a parabolic mass potential, 



± \ El 



+ sgn(?7 )V7 (2n+l) 



(92) 

The energy dependence with respect to the second dis- 
crete number n is now quite different from the previous 
Fock-Darwin-type spectrum in a scalar 2D-parabolic po- 
tential, Eq. ([86]) . As a specific illustration for the case of 
a circular parabolic mass potential V z (r) = (l/2)U~o(x 2 + 
y 2 ) together with a zero scalar term V s = 0, the discrete 
energy levels are clearly anharmonic with respect to n, 



E m ,n = -fU ± y/ (hfl c V^) 2 + [l 2 B U (m + n + 1/2)] 2 , 

(93) 

an expression which was already quoted in Eq. ([7|) . 



D. Discussion for arbitrary smooth potentials: A 
hierarchy of local energy scales 

It is worth emphasizing that for arbitrary two- 
dimensional potentials V s (R) and V z (R) that are smooth 
at the scale of the magnetic length Is, the present vortex 
formalism turns out to be extremely useful because it ex- 
plicitly puts forward the existence of a hierarchy of local 
energy scales. Such a hierarchy can then be exploited 
to devise successive approximation schemes, leading to 
controlled expressions for all physical observablcs at fi- 
nite temperature. This has already been proved with the 
concrete example of the temperature-broadened STS lo- 
cal density of states for the 2DEG (see Sec. IV of Rcf. 
[26l ). and the same mechanism holds also in the case of 
graphene studied here. 

To understand qualitatively the origin of this hierar- 
chy of local energy scales, it is useful to rewrite the *- 
bidiffcrcntial operator, Eq. (|51[) . under the equivalent 
form, 



1 / / 2 \ p - 
■ b \ (j V 



with 



p=0 y 



C = (d x ~$y - d Y d x 



(94) 



(95) 



The arbitrary large number of derivatives in expression 
([M| is clearly an indication of the nonlocal nature of 
quantum mechanics. However, and remarkably here, we 
realize that nonlocality manifests itself through quasilo- 
cality in the vortex representation. This is due to the 
fact that the nonlocal electronic Green's function G(r, r') 
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can entirely be determined from the knowledge of the lo- 
cal vortex function <? m (R), see connection formula (|54[) . 
This quasilocality property (which holds independently 
of the form of the potential landscape and thus can be 
seen as resulting uniquely from the coherent character of 
the vortex states), allows one to have a quasilocal quan- 
tization view. Clearly, the local Green's function g m (R) 
depends on the potential matrix elements v m (R) via the 
action of the * product, see Eqs. (|67|) - (|70l) . As obvious 
from expression (|M|) , each power of the bidifferential op- 
erator C acting on the functions u m (R) and g m (R) gen- 
erates higher and higher derivatives ^d^<? m (R) of the 
local Green's function associated with hierarchy of en- 
ergy scales of the type l^d^v m (R) ■ These energy scales 
get smaller and smaller at increasing p in the case of a 
potential smooth at the magnetic length scale, allowing 
one to control systematically the calculation. 

For instance, leading order expressions (jT5j) and ([75)1 
for the vortex Green's function were derived assuming 
that one can neglect potential curvature terms (associ- 
ated to the geometric invariants involving second-order 
spatial derivatives of the potential). This type of approx- 
imation is in fact controlled as long as temperature ex- 
ceeds the local energy scales appearing at next to leading 
order, respectively [7(R)] 1/2 of Eq. fl75J} and [?7(R)] 1/3 of 
Eq. (|79[) . In that case, quantum effects such as quanti- 
zation and tunneling are certainly missing, yet this ba- 
sic approximation already encodes the structure of the 
delocalized edge states far from the regions where the 
potential is strongly curved. 

We have seen in Sec. IIV CI that it is possible to go 
one step further by including the curvature contributions 
[term p = 2 in Eq. (|94j) ]. and this reintroduces quantiza- 
tion and tunneling in case of confined or open potentials, 
respectively. Again, one expects that the refined expres- 
sions obtained for the vortex Green's function [Eqs. ([76| . 
(|82j) . and ([88} depending on the dominant type of scat- 
terers] encode correctly the quantum dynamics down to 
further and even smaller energy scales associated to ge- 
ometrical invariants involving third order spatial deriva- 
tives of the potential. 

These considerations show the existence of a hierarchy 
of local energy scales formed by the successive spatial 
derivatives of the potential and hint that the passage 
from purely local physics (which is the hallmark of clas- 
sical mechanics) to highly nonlocal quantum-mechanical 
physics (which is the apanage of highly unstable quantum 
states) is worked out gradually when the temperature is 
progressively decreased. Therefore, at least in the large 
magnetic field regime, it is not needed to diagonalize nu- 
merically the random Schrodinger or Dirac equation in 
order to calculate precisely physical quantities, since tem- 
perature down to the Kelvin range in real experiments is 
not likely to be very small compared to the tiny energy 
scales at order l 2 B (for smooth potentials). What is ne- 
glected in our approximation scheme are contributions 
of some highly non-local quantum states superpositions, 
which are irrelevant in realistic experiments at finite tem- 



perature. 



V. CONNECTION WITH THE DEFORMATION 
QUANTIZATION THEORY 

A. Deformation quantization theory in classical 
phase space 

Before exploiting the expressions for the Green's func- 
tions derived in Sec. lIVl we would like to make important 
comments on the structure of the dynamical equations 
obeyed by the Green's functions g m (R) [general Dyson 
Eq. ([50)1 at any magnetic field or Eq. (|57|) in the absence 
of Landau-level mixing at high magnetic field]. After 
completion of Rcf. (2(| we have indeed realized that the 
★-bidifferential operator, Eq. (|5Tj) , involved in these lat- 
ter equations has a form analogous to the so-called star 
product, which has been the subject of intense research 
in mathematical and in high-energy physics because of its 
fundamental role in the principles themselves of quantum 
mechanics j^2r— More precisely, there have been many at- 
tempts to formulate quantum mechanics from a classical 
point of view, i.e., as a theory of functions on phase space, 
and one suggestion^ was to understand quantization as 
a deformation of the structure of (Poisson-Lie) algebra 
of classical observables. The ^.-deformation theory of the 
classical mechanics relies on the introduction of a star 
product, 



-kn = exp 



dp* - ^ f 



d x 



(96) 



in place of the usual product between phase-space func- 
tions. Here, x and p x are, respectively, the position 
and momentum which are canonically conjugate vari- 
ables. We discuss first here the quantization for a par- 
ticle in one dimension in the absence of a magnetic field 
(in two dimensions, classical phase space is four dimen- 
sional, see discussion in Sec. IVB[) . As a key princi- 
ple, the entire quantum dynamics is encapsulated in the 
noncommutative operator, Eq. (|96[) . which turns out 
to be the unique associative pscudodiffcrcntial deforma- 
tion of the ordinary product. Within the deformation 
quantization theory, the Poisson brackets of classical me- 
chanics between two phase-space functions f(x,p x ) and 
g(x,p x ) are replaced by the Moyal brackets^ defined 
as commutators (in the star-product sense) [/, g] m = 
(/ *h 9 — 9 *h f)/ifc- Obviously, Moyal brackets are H- 
dependent brackets which reduce smoothly to the Pois- 
son brackets in the limit h — s- (hence the origin of the 
"deformation" picture) . 

The deformation quantization approach appears as a 
generalization of original ideas put forward by Weyl, 
Wigner and Moyal^ (for a short historical account, see 
paper— and references therein), which were aimed at get- 
ting a sound insight into the correspondence principle 
between classical and quantum mechanics. The defor- 
mation quantization formulation has acquired a clearer 
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mathematical status 30 years ago with the work of Baycn 
et al. where its autonomous and alternative character 
with respect to other formulations of quantum mechan- 
ics, such as the conventional Hilbert space and path inte- 
gral formulations, has be en p roved (for the recent status 
of the theory, see Refs. [40 and |43j). Because the ba- 
sic continuous structure of the classical phase space is 
conceptually kept in the deformation quantization the- 
ory, classical mechanics is easily and transparently recov- 
ered via a smooth transformation, in full contrast to the 
conventional operatorial approach of quantum mechanics 
formulated in a Hilbert space (spanned by a countable ba- 
sis of square integrable states) where the emergence of a 
classical character from the quantum substrate appears 
singular and rather challenging. For this reason, it has 
been underlined^ that the deformation view is presum- 
ably the right way to look at quantization. 



B. Vortex Green's functions as a mixed 
phase-space formulation of quantum mechanics 

Now considering explicitly two-dimensional electronic 
quantum dynamics in the ordinary 2DEG, the standard 
deformation quantization theory introduces electronic co- 
ordinates {x,y) and momenta (p x ,p y ) as natural vari- 
ables in a four-dimensional phase space. In a large mag- 
netic field however, the electronic classical dynamics con- 
sists of a fast cyclotron motion, which is centered around 
a slowly moving guiding center R = [X, Y). In the pop- 
ular operatorial language of quantum mechanics, these 
two relevant degrees of freedom arc introduced by de- 
composing the electronic coordinate operator f = R + 17 
into a relative position 17 linked to cyclotron orbits and 
a guiding center position R = (X,Y). It is well known 
that the guiding center coordinate operators obey the 
commutation rule [X,Y] = il B , showing analogy with 
the canonical quantization rule between the position x 
and the conjugate momentum p x . Therefore, the square 
of the magnetic length, l 2 B , plays the role of an effective 
magnetic field-dependent Planck's constant. Moreover, 
cyclotron motion associated to the relative circular or- 
bits 17 leads to quantized Landau levels and at very large 
magnetic fields completely decouples from the guiding 
center dynamics. 

This physical discussion shows that the canonical de- 
scription of phase space in terms of electronic coordinates 
(x, y) and momenta (p x ,p y ) becomes awkward in a mag- 
netic field. Quantum mechanically, this is reflected by 
the property that states that are coherent with respect 
to both positions and moment a 44 ' 45 cannot be eigenstates 
of the kinetic part of the Hamiltonian associated to cy- 
clotron motion, contrary to the vortex states. With the 
benefit of hindsight, the program that we have followed 
in the string of recent papers^— is precisely the for- 
mulation of deformation quantization in a mixed phase 
space associated with the combination of discrete Lan- 
dau levels m and two-dimensional guiding center coordi- 



nates (X, Y), which correspond to physical space. For the 
2DEG, this decomposition is naturally encoded within 
the vortex states ^> m ^ of Eq. (jTUJ) , whose coherent 
character with respect to the guiding center R brings a 
doubly continuous parametrization of phase space, while 
the discrete quantum number is associated to a standard 
quantization of cyclotron motion. 

The general equation of motion at any magnetic field 
for graphene is then given by Eq. and simplifies into 
a dynamics in two-dimensional phase space (X, Y) given 
by Dyson Eq. (|67p in the large magnetic field regime, 
as cylotron motion giving rise to Landau levels exactly 
decouples from the vortex (or guiding center) motion. In 
that case, Dyson equation has precisely the form of a star 
product, see the obvious connection between the * oper- 
ator, Eq. (|5ip . of the vortex formalism and the *% prod- 
uct, Eq. (|96[) . of the deformation quantization theory. 
High magnetic field dynamics is thus isomorphous to a 
one-dimensional Schrodinger (for the ordinary 2DEG) or 
Dirac equation (for graphene) with conjugate variables 
X and Y. More specifically, if we consider the lowest 
Landau level (allowing one to forget the spinorial struc- 
ture proper to graphene) , Dyson Eq. (|6"T|) is equivalent to 
the standard operatorial formulation with the Hamilto- 
nian H = vq(X,Y), where the effective potential Vq(R) 
is given by Eq. (|G5|) . In that case, dynamics results from 
the commutation rule [-X",Y] = il B so that kinetic-like 
energy terms emerge from the identification of the conju- 
gate momentum to X with Px = HY/l B . We emphasize 
that this derivation is free of the ambiguities found in 
the path integral formulation^ and reproduces the low- 
est Landau projection method pioneered for the 2DEG 
by Girvin and Jachi 29 ' 30 The vortex formulation of phase 
space is however more general, because it allows to con- 
sider not only the projection onto arbitrary Landau levels 
at infinite magnetic field, but also the coupling between 
them for arbitrary magnetic field. 

Therefore, the scmicohercnt character of the vortex 
representation offers a local quantization view in high 
magnetic field because phase space reduces to the phys- 
ical space of guiding center coordinates R. When con- 
sidering the motion in complicated potential landscapes, 
this leads to the existence of a hierarchy of local energy 
scales, allowing one to describe smoothly the crossover 
from the semiclassical guiding center motion at high tem- 
perature to the fully quantum dynamics at very low tem- 
perature, as discussed in Sec. IIVD1 



VI. LOCAL DENSITY OF STATES 
A. Generalities 

We now use the formalism developed in the previous 
sections and the resulting expressions for the graphene 
Green's function to investigate the characteristic features 
of the local density of states (LDoS). The goal of Sec. 
IVII is to show that a lot of information concerning the 
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different potentiais at piay in graphene can be extracted 
from the widths and shapes of the LDoS peaks in a high 
magnetic field. 

The LDoS is related to the electronic Green's function 
via the formula 



p(r,oS) = ImTrG(r, r, oj). 

TT 



(97) 



Note that with Eq. ([54]) . we can directly write the 
LDoS in terms of the modified local Green's func- 
tion 5mi,Ai;m2,A 2 (R-)- I n the case where the modi- 
fled Green's function is diagonal with respect to the 
Landau-level quantum number, i.e., g mi ,\iwi 2 ,\2 (&) = 
(R), we have to evaluate the action of 
the exponential differential operator onto the product of 
two vortex functions with identical Landau level, as done 
in Eq. ([63]). We therefore find that the LDoS [Eq. (|97| ] 
can quite generally be written in the absence of Landau 
level mixing as 
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where the kernel A rn (R) has been previously obtained 
in Eq. ([63]). We have also taken into account here the 
spin and valley degeneracies, which provide an overall 
prefactor of 4 when evaluating the trace in formula (1UT1) . 

In actual experimental conditions, one never has a di- 
rect access to the zero-temperature LDoS, due to an ex- 
trinsic smearing occasioned by the finite temperature T. 
The STS spectra at fixed energy e are proportional to 
the temperature broadened LDoS 
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formulas (|98|) and (|99|) . we get 
rd 2 R 



n' F (£-£ (R + r))A (R) 



+ 2 E E ~ W( R + r )) {C 1 + e ^»( R + r )) 

m— 1 e— ± 

x^ m (R) + (l-6/3 m (R + r))A m _ 1 (R)} ], (100) 

where the effective energy £ m ,±(R), the electron-hole 
asymmetry parameter /3 m (R) and the kernel K m (R) are 
given respectively by Eqs. ([74"]) . (jM)) and (jSB")) . The ker- 
nels K m (R) are oscillating yet normalized functions that 
are localized around R = on a characteristic length 
scale L m = IsV^m + 1, which one associates with the 
cyclotron radius. Only for the lowest Landau level m = 
does this length reduce to the magnetic length Is- 

In principle, one cannot strictly set the temperature to 
zero in Eq. (|100[) unless the effective potentials u^(R) 
and u~ (R) which compose the function £ mjE (R) are glob- 
ally flat. Indeed, for arbitrary potentials u+(R) and 
u~(R), it is important to have in mind that expression 
(|100P overlooks the fine structure of the zero-temperature 
local density of states, which requires to take into ac- 
count all existing spatial derivatives of these potentials 
^to(R) [see Eq. (|57)) ]. Nevertheless, it captures accu- 
rately the shape of the LDoS when the temperature ex- 
ceeds the (smaller) energy scales involving second and 
higher derivatives (in orthogonal directions) of the po- 
tentials w^(R) associated to curvature. Basically, under 



the inequalities L m \ V£ m , c (r)| > T ^> y |7m(R)| ; one ex- 
pects that the temperature gives a small contribution to 
the smearing of the LDoS in comparison to the intrinsic 
smearing generated by the spatial dispersion of the func- 
tion £ mj£ (R + r), (i.e., by the potential gradients) when 
performing the integration over the vortex position R in 
Eq. (1TTKJ1) . 



2. High-temperature regime 



where n' F (oj) = — l/[4Tcosh (u;/2T)] is the derivative of 
the Fermi-Dirac function. 



B. LDoS for locally flat potentials 

1. General expression 

The leading order result for the vortex Green's func- 
tion, Eqs. (fT3"| and (|75|) . applies when the disorder po- 
tential is locally flat on the scale I b ■ Mathematically, this 
approximation is controlled for temperatures larger than 
the smaller energy scales associated to local Gaussian 
curvature, such as Eq. (|78[) . In that case, using previous 



At very high temperatures such that T 3> 
L m \ V£ TO)e (r)|, the spatial dependence on the vortex po- 
sition R inside the Fermi derivative function can be ne- 
glected [here, we also disregard the R dependence of the- 
smooth function /? m (R + r)], so that expression (jlOOp 
simplifies into: 
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This semiclassical expression provides LDoS peaks of 
width 2T that are centered around the effective Landau- 
level energies £ m ,±(r) given by Eq. (f?4")) . In this regime, 
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the thermal broadening of the LDoS peaks is thus in- 
dependent of the Landau-level index, and the electron 
and hole peaks are characterized by the same heights. 
At lower temperatures, we now show that different 
linewidths, line shapes and particle-hole asymmetries are 
generated in the LDoS spectra, providing additional in- 
sight into the underlying scalar and mass potentials. 



first assume here that the potential is well approximated 
by its first-order gradient on the whole cyclotron orbit 
of radius L m , i.e., £ m , e (R + r) ~ £ m ,e(r) + R.V£ m , e (r). 
We can then perform analytically the Gaussian integral 
over R in Eq. (|100[) and obtain the intuitive result for 
the zero-temperature LDoS (see Appendix [D]) , 



3. Low-temperature regime for potentials smooth on the 
cyclotron radius L m 

In case when T < L m |V£ m)e (r)|, the spatial depen- 
dence of the Fermi function derivative must be kept. We 
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with F^ e (r) = /s|V^ m!e (r)| the local energy scale as- 
sociated to the drift motion and H m {x) the mth Hcr- 
mite polynomial. In order to keep the above expression 
compact, we have only written the zero-temperature lo- 
cal density of states, but the STS local density of states 
is readily obtained from Eq. (|99[) . The above expres- 
sion is quite reminiscent of the expression that can be 
obtained with the usual Landau states, of course gener- 
alized to the two-component spinorial structure proper 
to graphene, and taking into account that the potential 
landscape varies slowly in space [obvious from the r de- 
pendence of the width T^ e (r)]. In Sec. PvTC! we will 
further analyze expression (|102j) when discussing recent 
STS experiments. 

We note yet that for a given disordered potential land- 
scape Eq. (|102[) breaks down for sufficiently large quan- 
tum numbers m, because very wide cyclotron orbits of 
radius L m rj \f2mlB 3> Ib may explore random spa- 
tial variations in the potential. In that case, the more 
general expression (jlOOp is still valid, provided that the 
potential is smooth on the smaller scale I b (this is always 
the case at high enough magnetic field). This regime is 
now investigated. 



Low-temperature regime for potentials with random 
spatial fluctuations on the cyclotron radius L m 



spatial variations along the trajectory remain however 
smooth at the smaller scale Ib, general expression (|100[) 
for locally flat potentials is the one to consider. In order 
to get some analytical insight, we compute here a disor- 
der averaging of the LDoS. This procedure is clearly valid 
in two cases: (i) for the LDoS at very large Landau in- 
dex m 3> 1 , as very wide cyclotron radius L m can explore 
random configurations of the scalar disordered potential 
V s (r). Because of the large quantum numbers involved 
here, one should recover a semiclassical limit, as we will 
see; (ii) for any m and finite magnetic length (the fully 
quantum regime) , if one rather considers the sample aver- 
aged density of states (DoS). We stress beforehand that 
the LDoS at small m does not show self-averaging. In 
both situations, the computed averaged density of states 
is a spatial-independent quantity. The calculation per- 
formed in Appendix iDl provides the following result: 
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with the characteristic energy width L^ given by 



In cases where the disorder potential fluctuates spa- 
tially on the scale of cyclotron radius L m , formula (|102l) is 
clearly invalid, as the linearization of the effective vortex 
potential £ mi£ (R + r) cannot be made anymore. When 



(104) 
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where S(q) is the Fourier transform of the potential corre- 
lation function (see Appendix iDl for details) and K m (R) 
was denned in Eq. (|63[) (we write here K—i = Kq in or- 
der for the above formula to apply at in = as well). In 
order to simplify the derivation, we have assumed that 
the antisymmetric part V z of the total potential V can 
be neglected compared to the diagonal scalar component 
V a . 

Equation (| 104|) can be first analyzed in the following 
semi-classical limit, Ib —> and m>l, while keeping the 
cyclotron radius L rn = y2m + 11 b fixed. In that case, 
the function K m (r), which is peaked at the distance |r| = 
\J2vhIb — L m with a width Ib, becomes a delta function 
along the cyclotron radius, K m (r) ~ 2 ^ L 5{\r\ — L m ). 
In this semi-classical regime, we recover results derived 
by other means 2 ^ for the 2DEG, namely, 

[T » oS? = f ^ 2S{q)Vh[Lmq)? ^ M^M x i t 

J 2tt n-L m y/m 

(105) 

where the asymptotic limit of the zeroth order Bcsscl 
function Jo was used, assuming the disorder to be ran- 
dom on the scale L m , so that the integral in Eq. (|105[) 
is dominated by its tail. We note that our expres- 
sion (|104j) is more general than the above result (|105l) . 
because it also describes the averaged density of states 
for any m (including the strong quantum regime at finite 
Ib). Clearly, our calculation incorporates wave function 
spreads on the scale Ib, a purely quantum length scale 
which has completely disappeared from the semiclassical 
result [Eq. (|105[) ]. In all cases (semiclassical or quantum 
dynamics), the general trend is that the cyclotron motion 
averages out the local potential at increasing radius L m , 
so that the width of the DoS decreases with m. This ef- 
fect is discussed now in more detail at the light of recent 
LDoS measurements. 



C. Interpretation of the STS experiments 

Recent experimental works by Li et alJ^ and Miller 
et alw^ have investigated by STS the LDoS in graphene 
at high magnetic field and have revealed the relativistic 
nature of the Landau levels in the measured energy spec- 
trum. Besides this precise verification of the sequence of 
graphene Landau levels at the energies +Ml c y/m, one can 
note several other striking aspects of the data. At a given 
large magnetic field and for a fixed tip position, the width 
of the mth Landau-level peak in the LDoS p STS (e,r) 
is seen to grow as \fm with increasing Landau-level in- 
dex m, as demonstrated in Ref . [TT| and also observed in 
Ref. [l2|. At the same time, the LDoS peaks display an 
energy dispersion MS cl function of tip position, reflect- 
ing the underlying effective potential, see the discussion 
in Ref. [TH Quite contrary to the fixed tip LDoS peaks, 
the energy spread of the spatially averaged mth Landau 
level decreases with m. This effect is easily understood 
on general grounds by the smearing of the local poten- 



tial by larger and larger cyclotron orbits, as discussed 
above and embodied in the DoS width r^° 5 of Eq. (|104|) . 
In particular, the semiclassical limit (m ^> 1, Ib — > 0), 
which does not completely apply to the experiment for 
which Landau levels are only observed up to m = 7, gives 
the result T®° s cx to -1 / 4 , as first derived by Raikh and 
Shahbazyar*2£ for the non-relativistic 2DEG, showing a 
clear decrease in the width with m. 

As an illustration of truly quantum smearing of the cy- 
clotron motion at finite Ib for the first few Landau levels, 
which corresponds more to the actual experimental situ- 
ation at high magnetic fields, we have plotted in Fig. Q] 
for m < 4 the effective potential in graphene obtained 
from Eqs. and l[74"|) in the case of negligible band 
mixing [i.e., |v~(r)| < M2 C ], 

^ m>+ {r)=E m + - J d 2 r l V s {r 1 )[K m {v-r 1 )+K m _ 1 {r-r 1 )] 

(106) 

as a function of tip position r and for a given (scalar) dis- 
order realization, obtained as a superposition of localized 
long-range potentials. The upper panel of Fig. [TJ which 
corresponds to a (uni-dimensional) disordered scalar po- 
tential landscape V s (r) smooth on the scale Ib, shows 
that the effective potential £ mi +( r ) follows precisely the 
bare disorder potential for the lowest Landau level m = 
0, yet presents some moderate deviations for the fol- 
lowing levels, illustrating the small averaging present 
on the larger scale of the quantum cyclotron radius 
L m = \J1m + 11b- In contrast, the lower panel of Fig. Q] 
presents the situation of a disordered scalar potential 
landscape V s (r) which has spatial variations comparable 
to Ib [we stress again that the effective potential given 
by Eq. ([73|l and thus also by Eq. (|106[) has a truly non- 
perturbative character in lg]. In that case, we can no- 
tice two effects: (i) the effective potential £ m-+ (r) shows 
important quantitative deviations from the bare one al- 
ready in the lowest Landau level m = 0; (ii) at increas- 
ing m > 0, stronger and stronger averaging effects take 
place, so that the effective potential £ mi +(r) rapidly flat- 
tens out. As a consequence, the typical energy width of 
the effective potentials £ TO ,±(r) as a function of position 
r clearly decreases with growing m. This effect is clearly 
seen in the STS data of Ref. [jjg for graphene and can 
be also recognized in recent measurements on standard 
2DEGs by Hashimoto et alZ 

We now discuss in more detail the STS spectra taken 
at fixed tip position, presented in the experimental pa- 
per a n i 12 that showed a broadening of the Landau levels 
with a ypm scaling at increasing m. At high temper- 
atures, such that T ^> L m | V£ m ,±(r)|, the broadening 
has a purely thermal origin, with a fixed width set by 
T and an exponential line shape (given by the Fermi 
function derivative). It is worth noting that the appar- 
ent increase with m of the heights of the LDoS peaks in 
graphen o 11 ' 12 is solely due to the collapse of Landau lev- 
els, E m+ i — E m oc 1/y/rh at large m, yet the underlying 
Landau peaks show a width insensitive to m. 

Contrary to the discussion given in Ref. [Til we empha- 
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FIG. 1: (Color online) Dimensionless effective potential 
£ mj +(r)/M7 c = y/m + v^i(r)/Ml c from Eq. (|106|l as a func- 
tion of linear tip position t/Ib for the first Landau levels 
m = 0,1,2,3 (bottom to top in full lines), and compared 
to the bare potential energy sjm + V s (r)/hfl c given by the 
dashed lines. The top panel corresponds to smooth disorder 
while the bottom one has stronger variations in the potential 
on the scale Ib (see the relative axes). 



size that results of disorder averaged density of states, 
such as our Eq. (|105|) or the formula obtained, e.g., 
in Ref. H^, do not apply to account for the width of 
the STS peaks at fixed position, for which an expres- 
sion for the local density of states, such as Eq. (jlOOp 
or (|102[) . should instead be considered. In fact, the en- 
ergy spread of the Landau-level peaks observed exper- 
imentally at low temperature in the LDoS can be eas- 
ily understood to originate from wave-function broaden- 
ing. Indeed, in formula (|102[) for instance, the polynomial 
\H m (x)\ 2 being of order 2m, the squared wave function 
f(x) = \H m (x)\ 2 e~ x turns out to be spread on a char- 
acteristic scale x m = \j2m + 1. We note that in Fig. 2 
of Ref. H(| a square-root dependence of the Landau-level 
widths with the Landau-level index can also be observed 
at zero-temperature (the oscillatory substructure of each 
Landau level peak disappears when including a small 
thermal smearing, as performed here). Turning to the 
LDoS expression (|102[) . one sees that the effective energy 
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FIG. 2: (Color online) Energy-dependent STS spectra for 
the local density of states p STS (r, e,T) at fixed tip posi- 
tion r from Eqs. ([99]) and (|102|) in units of 4/(27r/|) and 
as a function of energy e for several temperatures T/Kl c = 
0.03,0.05,0.07,0.09 (top to bottom). Here the local energy 
scale associated to the drift motion in the lowest Landau level 
is To oc (r) = l B \W s \ = 0.02fif2 o . At the lowest given temper- 
ature, the smearing with increasing m of the Landau-level 
peaks is characterized by a local energy width which roughly 
grows as \fm (thermal smearing provides still some dominant 
contribution). 



width of the mth Landau-level peak is roughly given by 
the local energy scale \/2m + 11b\ VV^(r)| = L m |VV^(r)| 
[here we have used the fact that the effective poten- 
tial £m,±( r ) roughly follows the bare potential V^(r)], 
which scales as sfm as observed in the experiment ^ 
This effect can be checked by a simple numerical eval- 
uation of Eq. (|102p , taking into account the convolution 
with a thermal smearing as resulting from Eq. for 
the STS local density of states. The obtained result 
for the sequence of LDoS peaks is shown in Fig. [2] for 
different temperatures. At temperatures comparable to 
Fg oc (r) ~ is|VV^|, quantum smearing due to the drift 
motion, which is encoded by the spatial dependence of 
the kernel K m (R) in the general expression (|100|) . or by 
the Hermite wave functions in the special case of a glob- 
ally flat potential [see Eq. (|102[) ]. starts to appear. The 
growth of the energy width of the LDoS peaks at increas- 
ing Landau-level index m is visible for the lowest chosen 
temperature in Fig. [2J In that case, one also sees a clear 
decrease in the heights of the LDoS peak with m, as ob- 
served experimentally ^ 1 1 1 12 Because the total smearing 
of the Landau levels depends both on thermal and in- 
trinsic wave function broadening, the linewidth is only 
roughly behaving as \]rn. 

Finally, we address the question of the Landau-levels 
line shape in the LDoS. In experiment of Ref. [ill, it has 
been pointed out that Lorentzian fits are significantly 
better than Gaussian ones to account quantitatively for 
the broadening of the LDoS peaks. On the other hand, 
in experiment by Miller et al*£ the line shape has been 
modeled by a convolution of Lorentzians and Gaussians 
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to include extrinsic origins of broadening induced by 
temperature and instrumental resolution. On theoreti- 
cal grounds, thermal broadening implies exponential line 
shapes (in between Lorentzians and Gaussians), while in- 
trinsic wave function broadening of drift states (for non- 
vanishing local potential gradients) leads to Gaussian- 
type energy decay. We also note that spectra taken in 
regions of small potential gradients involve intrinsic ex- 
ponential linewidth due to curvature effects, see Ref. [3l| 
for a discussion of the lowest Landau-level LDoS peak in 
the 2DEG. Therefore, it is difficult in general to disen- 
tangle the different contributions from the experiment, 
and systematic studies in temperature and as function 
of tip position, would be required to settle precisely this 
issue. 



VII. CONCLUSION 

In this paper, we have extended to the graphene case 
a Green's-function formalism well suited to study the 
mechanism of lifting of the Landau-level degeneracy by 
a smooth potential landscape at high magnetic fields, 
which was originally developed for the two-dimensional 
electron gases. The whole formalism relies on the use 
of a particular representation of semicoherent states, 
which are eigenstates of the kinetic part of the Hamil- 
tonian. These so-called vortex states in the 2DEG case, 
or graphene vortex states in the graphene case, are both 
characterized by an integer topological quantum number 
to, related to the vortex circulation and giving rise to 
the Landau quantization of the orbital motion, and by 
a doubly continuous quantum number R, corresponding 
to the location of the vortexlike phase singularities of 
the electronic wave function and characterizing the huge 
degeneracy of the Landau levels in the absence of disor- 
der. The coherent states character with respect to the 
degeneracy quantum number R allows one to project the 
electron dynamics onto this overcomplcte representation 
of states, which rigorously extends to quantum mechanics 
the classical guiding-center picture. 

In a first stage, we have derived the exact matrix el- 
ements for smooth arbitrary scalar and mass potentials, 
as well as for off-diagonal smooth potentials related to 
ripples in graphene. The particular form of these matrix 
elements has revealed the different processes leading to 
Landau-level mixing and coupling between electron and 
hole bands. We have shown that at high magnetic field, 
when the Landau-level mixing can be safely neglected, 
a mixing between the hole and electron energy bands is 
unavoidably induced by second-order derivatives of the 
scalar potential, independently of the presence or not of 
a mass potential. We have been able to derive in this high 
magnetic field regime exact expressions for the electronic 
Green's function in the presence either of an arbitrary 
quadratic scalar potential or an arbitrary quadratic mass 
potential. 

Besides affording the derivation of unique Green's 



function solutions valid for closed and open quadratic 
potentials which underline the dual correspondence be- 
tween quantization effects and tunneling effects, we have 
emphasized that the semicoherent vortex representation 
offers a quasilocal perspective of the quantization process 
closely related to the deformation view of the classical 
phase-space mechanics, a property which turns out to 
be essential to capture the transition from the nonlocal 
quantum world to the local classical world. Furthermore, 
the vortex representation has revealed a hierarchy of lo- 
cal energy scales formed by the successive derivatives of 
the potential and thus ordered by their degree of nonlo- 
cality. As a result, quantum features associated with the 
lowest derivatives of the potential appear to be the most 
robust against the inelastic effects. We have emphasized 
that the consideration of a finite temperature allows one 
to disregard the smallest inaccessible energy scales and 
thus to devise successive approximation schemes for an 
arbitrary smooth potential. 

Within this spirit, we have derived controlled analyti- 
cal expressions for the local density of states in graphene 
valid at high magnetic field in the presence of smooth ar- 
bitrary scalar and mass potentials within different tem- 
perature regimes. We have identified the most relevant 
mechanism of intrinsic broadening of the LDoS peaks 
and have shown that a lot of information on the dif- 
ferent potentials at play in graphene can be extracted 
from the experimental LDoS spectra performed at high 
magnetic field. Finally, we have been able to explain 
a few of the experimental findings, e.g., concerning the 
scaling of the LDoS peaks with the Landau-level index, 
recently observe d 11 ' 12 in scanning tunneling spectroscopy 
of graphene. 
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Appendix A: Energy spectrum for closed quadratic 
potentials 



In this appendix, we show how the energy spectrum for 
a confining quadratic potential (with a positive Gaussian 
curvature) can be determined from a retarded Green's 
function expression such as given, e.g., by Eq. (|76[) . For 
7o(R) > 0, the function 



W(R,t) 
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is periodic in time with the period T — 2ir/ (R) at 
fixed R. We thus expand it in a Fourier series 



W(R,t)= a p (R)e- ip V^( R ) t , 



(A2) 



and insert expression (| A2|) into Eq. ([76]) to straightfor- 
wardly get after integration over time 
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with tco(R) = 5o(R) - 77o(R)/ 7o (R). 

The Fourier coefficients a p (R) are given by 
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(A5) 



with p(R) = t/o(R)/ [7o(R)] 3 ^ 2 - We rewrite the following 
function appearing in the integrand of integral (|A5|) as 



meaning that the series on the right-hand side of Eq. 
(|A8j) is divergent for x < and z — e txp . 

Using Eqs. ([A6]) or (jA7| . and Eq. (jA8|) by writing 



2|p(R)| and z 



depending on the sign of 



the quantity p(R) [we introduce the short-hand notation 
X = sgnp(R)], we can easily perform the integration over 
the angle 9 in Eq. (|A5|) and find 



Op(R) = 2(-l)"e-W R )li„ (2 |p(R)|) (All) 



for p = x(2n + 1) and a p (R) = for any values of p ^ 
X(2n + 1). Therefore, only the terms with p = x(2n + 1) 
remain in expression (|A3j) . where n is a positive integer 
and x = ±1 is an index determining if the region is locally 
convex or concave. 

Now, for the particular case of purely quadratic scalar 
and mass potentials, the poles of the Green's function 
(|A3[) arc R independent, and thus directly yield the en- 
ergy spectrum, Eq. (|81[) . with the set of quantum num- 
bers (m, n) if the quadratic potential V s — V z is convex, 
i.e., confining (x = +1 in this case). 
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Appendix B: Solution for a locally quadratic scalar 
potential V s 



It is then convenient to introduce the identity (see, e.g., 
Ref. Et" 



1 / xz \ + °° 



(A8) 



n=0 



where L n (x) is the Lagucrrc polynomial of degree n. For- 
mula (|A8p is usually defined for |z| < 1, but it can be 
checked that it still holds for z = e %v with ip ^ 2irj (j 
a positive or negative integer) at i > 0. Indeed, using 
the asymptotic behavior of the Laguerre polynomials at 
large n and x > 0, 
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we note that the series on the right-hand side of Eq. (|A8|) 
is semi-convergent (this can be established using Abel's 
test). On the other hand, for x < 0, we have 



L n (x) 




20F(n|x|) 1 / 



In this appendix, we solve the equations of motion, 
Eq. ([67]) , in the regime where we can consider that the 
effective potential u~(R) has a negligible spatial disper- 
sion and that the effective potential £>+(R) can be locally 
described up to its second-order derivatives (i.e., it is lo- 
cally written as a two-dimensional quadratic potential). 
These assumptions turn out to be exactly fulfilled in the 
particular case of a globally quadratic scalar potential 
V s (R) and a constant mass term V z , for which the ma- 
trix elements at high magnetic field read for m > 1, 



5+(R) = V s (R) + m l -fA n V s (R) 



(Bl) 



v~(R) = -Vg + -fA R V s (R) = est. (B2) 



Using the explicit form, Eq. (15TT) . of the * operator and 
the fact that ?)+(R) is a quadratic function (so that all 
its derivatives higher than 3 vanish) and that £>~(R) is 
quasi-independent of R, Eq. (|67[) becomes 
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I 2 

(bJ - U^j(R) - Em,*! + «0 + ) <?ro;Ai;A 2 (R) = #Ai,A 2 + V5m;-A i; A 2 (R) + «y [cbfU+dy - SyU+dx] 3m;A i; A 2 (R) 

-f 9 2 x + (d 2 x i+) d 2 Y - 2 (d x d Y v+) dxd Y ] S™ ; a i; a 2 (R). (B3) 

To solve Eq. (|B3|) . wc introduce an arbitrary reference point Ro and write <7m;Ai;A 2 (R) = fm;\ 1 ;X 2 [E(R-)] with 
E(R) = u+(R) — u+(Rq). Substituting this form into Eq. (|B3[) . we get the system of differential equations, 



(Tm^ + ^™)^2 + - E + CJ - u+(R ) - + i0+ 



fm;Xi;X 2 (E) ~ V m f m --\ i; \ 2 (E) — #Ai,A 2 ; 



(B4) 



where the geometric coefficients 7+ and 77+ have the same definitions as in Eqs. (|78|) and ([79| [with vo(R) replaced 
by the effective potential v+ (R)], and are expressed at the reference point Ro. In the derivation of Eq. (|B4p . we have 
used the relation ry„(R) = T)+ + 7+ [u+(R) — u+(Ro)] which holds for any quadratic potential S+(R). We then go 
to Fourier space by writing 



fm;Xi;X 2 ( E ) = / dTF m-M-M(j) 



(B5) 



and obtain a system of coupled first-order differential equations for F 



^to;Ai;A 2 (t) 
VF m; _A i; A 2 (T) = 5(t)S Xu x 2 - 



(B6) 



Introducing into Eqs. (|B6[) the (last) change in function, 



F 7 n - , 



/ x _ ^Ai;A 2 [^( r )] i[«-e+(Ro)]t(T) 
i|Ajv J r 1 — r e 



7m T 



with 



t(r) 



xe 



arctan 



i(^/7+)[t(r)-T] 



(B7) 



7m ^ 



(B8) 



we arrive at a simple linear system of two coupled first- 
order inhomogeneous differential equations with constant 
coefficients, 



dt 

; d 

'dt 



with 



XE„ 



\E n 



70 H 



70 H 



of Eqs. (|B9[) and (|B10|) leading to a well-defined integral 
(|B5|) can then be readily derived and reads 



h x -x(t) 
hx-,-x(t) 



i9(t) 
2 

f — Xa 

-0m 



I + \a n 

0m 



- it V / - E ™ + [ ii " 



f .(BI2) 



hx ; x(t)-v-h-x-x(t) = S[r(t)}, 



h-x-,x(t) — v m hx ; x(t) = 



The expressions for the coefficients a m and j3 m are given 
in Eqs. and Coming back to the original func- 

tions <7m;Ai;A 2 (R) an d setting R = Ro, we get the com- 
pact expression for the modified retarded Green's func- 
tion written in Eq. (|82|) . For the case 7+ > 0, it is im- 
portant to realize that the relevant variable is the time 
t, not the variable r (whereas it is possible to work in- 
(B9) differently with t or t for 7+ < 0). It can be checked 
that integral (|82p is well defined as well for 7+ < as 
for 7+ > [in the latter case the infinitesimal quantity 
' i0 + is crucial while it does not help to make the integral 
convergent when expressing the solution under the form 
of an integral over r as within Eq. (|B5p ]. 



r(t) 



tan 



(Bll) 



Note that S[r(t)] = S(t) if 7+ < and 5[r(t)] = 
J2 n S (t — mt/x/^ij if 7+ > 0. Let us consider for the 
time being the case 7+ < 0. The solution of the system 



Appendix C: Solution for a locally quadratic mass 
term V z 



In this appendix, we solve the equations of motion, 
Eq. (|67p . in the regime where we can consider that the 
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effective potentiai u+(R) has a negligible spatial disper- In the regime considered in this appendix, Eq. ([67]) bc- 

sion and that the effective potential v~(R) can be lo- comes 

cally described up to its second-order derivatives. This 

regime contains as a particular case the situation where 

the scalar potential V S (R) is globally constant in space 

and the mass potential V z (R) has a quadratic dependence 

on R. In this particular case, we obviously get exactly 

that ?)+(R) = est and u~(R) depends quadratically on 

the variable R, 

S+(R) = V S (R) - ^A R V Z (R) = est, (CI) 
S-(R) = -V z (R)-m l fA R V z (R). (C2) 



E mM + i0+] 



In 

<? m; A i; A 2 (R) = 5 Xl ,\ 2 + Wm( R )ffm;-Ai;A 2 (R) + «y [dxV~d Y ~ d Y V~d X ] g m; -A i; A 2 (R) 

-f [( d YV~) d\ + (d 2 x v~) d 2 Y - 2 (9 X ^«-) 9x«9y] ff m; _ AliA3 (R).(C3) 
I 



As in Appendix IB1 we introduce a reference point Rq. 
It can then be guessed that the functions g m ;Ai;A 2 ( R ) arc 
functionals of the potential -E'(R) = v r ~ (Ro)— v™(R), i.e., 

we can write <?m;Ai;A 2 ( R ) = fm-MM [ E Qty]- Tnc con- 
tributions in Eq. (|C3[) involving the first-order deriva- 
tives of the function w~(R) then vanish. Furthermore, 
we shall suppose that the equality (f7T|) still holds in the 
present studied case, what can be justified a posteriori. 
The problem thus reduces to the resolution of a system 
of two coupled differential equations. Indeed, Eq. (|C3[) 
yields the system of equations 



with 



d 2 d _ , 



fm;-\ 1 ;\ 2 {E) 

-{lo- v+ - E mM + iO + )f m;Xv ,x 2 (E) = -S XlM (C4) 



s(r) 



arctan 



7m T 



(C7) 



with the coefficients 7" and 77" expressed at the posi- 
tion Ro and given by the formulas (|78|) and (|79|) written 
for the potential — u m (R) in place of vq(R). Applying 
the Fourier transformation, Eq. (|B5|) . we arrive at the 
system, 



a simpler system of two coupled differential equations 
with constant coefficients comes out, 



i—h-\ ; \(s) - (ui - u+ - XE m + iO + )h x -\{s) 

= S[r(s)}, (C8) 

i-^-h x . x (s) -{lo- v+ + \E m + iQ+)h_ x , x {s) = 0. (C9) 
as 

Note that, in contrast to the situation encountered in 
Appendix [Bj the variable s does not have the meaning 
of the time here since it is no more conjugated to the 
frequency lo [this is the reason why we took care of nam- 
ing the variable differently here although the expressions 
([Blip and (|C7[) are almost identical]. After diagonaliza- 
tion of the 2x2 system, we obtain that the homogeneous 
solution of Eqs. ([U8j and (JUS]) is 



K 1 + ^™ T2 )^ + *T™ T ~ VrrJ 2 + £ TO (Ro) 



F, 



m: — \ 



i;A 2 (t) 



-{lo - 5+ - E mM + iO + )F m . Xl .x 2 {T) = -5{t) 8\ lt \ a 



h\-\{s) 
h-x ; x{s) 



C 



-D 



(C5) 



+ \E n 
lo-v++ XE m 

Km 



(CIO) 



Introducing into Eq. (|C5|) the change in function 



F n 



,s = h\i;\ 2 [s{t)} p i[v-(Ro)+r,-/-y-}s(T) 



\/l + 7m T 



xe 



with C and £> two arbitrary constants, and the energy n m 
given by Eqs. (|59")l and ([9"0"|) (we forget for the time be- 
ing the infinitesimal quantity i0 + ). The inhomogeneous 
solution of system of Eqs. (|C8[) and (|C9[) is then ob- 
tained by varying the constants C{s) and D{s). As a 
i(flm/7^)T ( C6 ) result) we get C /( S ) = = -i<5(s)/(2/t OT ), that is 
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C(s) = Ti9(±s)/(2K m ). Using that 

(C 

the sign ± for the functions C(s) and -D(s) is then chosen 
in such a way that integral (|C11|) is convergent with the 
help of the infinitesimal quantity i0 + . Finally, taking the 
reference point Ro = R [so that E(R) = 0], we arrive at 
the expression (|88p for the Green's function, which holds 
irrespective of the sign of the coefficient 7~(R). 

Appendix D: Simplifying the LDoS expression for 
locally flat potentials 

In this appendix, we simplify further the expression 
(IIOOP for the LDoS (valid for locally flat potentials) in 
the low-temperature regime within two different cases: 



(i) case of a potential landscape which varies slowly on 
the scale L m (Appendix lD (ii) case of a potential land- 
scape which fluctuates spatially in a random way on the 
scale L m ( Appendix ID 21) . 



1. Potentials flat on the scale L m 

Writing the derivative of the Fermi-Dirac function as 

(Dl) 

and using the linearization of the effective energy 
£m,e(R + r) ~ U,e(r) + R.V£ m , e (r) in Eq. (HOOj) , we 
can then perform the Gaussian integral over the vortex 
position R to get the LDoS expression, 



p^(r, £ ,T) 



1 4 



2ttZ| tt 



+ 00 

dt 

-OO 



e *t(e-«o(r)) exp 



+ (1 - ep m ) 



i d 



m—1 



t 2 rf, oc (r) 2 
i 



i + °° r 

EE 

Z m=le=± L 



1 d r > 
nl ds 



(m - 1)! ds" 1 - 1 J 1 - s 



e^-WM) exp ■ 



i 2 r^ e (r) 2 



4A, 



(D2) 



with r^ c e (r) = Ib\ V£ mj£ (r)| and assuming temperature is low enough [i.e., T <g; r^ c e (r)] so that the limit T — > can 
be taken. We then perform the integral over time t, and obtain: 



1 4 



1 



T loc {r) 

+(1 - ep m ) 



exp 



r' oc (r) 

(m - 1)! ds" 1 ' 1 J 7T 



, } + 5££i5w{ (1+ ^ 



1 d" 1 
n! <9s^ 



: exp < -A s 



£ - £m,e(r) 
pioc C r -j 



s=0J 



(D3) 



Finally, using the following relational obeyed by the Her- 
mite polynomials H n {x) 

V L 1 J „_Q 

formula (|D3j) can be recast in expression (|102|) . 

2. Potentials random on the scale L m 

We consider here the limit where the potential has 
strong spatial variations along the cyclotron radius, 
which applies to the situation of large Landau levels. We 
assume for simplicity that the antisymmetric part V z of 
the disorder potential can be neglected compared to the 



scalar component V s so that the effective potential given 
by Eq. (J73) reads 

f m ,±(r) = ±E m + \ J d 2 V V s ( V )[K m (R-r 1 ) 

+A m _i(R-7 7 )], (D5) 

where K—i = Kq in order for the above expression to 
apply for m = as well. 

The averaging procedure is carried through the 
isotropic distribution function S(q) in Fourier space (here 
q = |q|) that describes the spatial correlations of disorder 

(D6) 



2G 



so that the spatially averaged LDoS becomes 

a r+oo /vj2t> +°° -i 

x y VV S [K m (R - r) + *T m _i(R - r )] e ^-?^W]* 

xexp|-i y"d 2 R! Jd 2 K 2 S-^Ri-RaJV.CEliJViCRa)! 

(D7) 

where the distribution obeys <5(R) = J" d 2 r; <S _1 (R — 
77)5(77). Inserting the effective potential, Eq. (|D5|i . and 
performing the functional integral over the disorder real- 



izations, we obtain 

17 B rn=0 e=± 

x [tf m (R - r) + tf m _i(R - r)] cxp j-Jt 2 [r£ oS ] 2 J 

(D8) 

where the width r^ oS is given by Eq. (|104[) . The above 
expression has obviously become r independent so that 
the R integral can be carried using the normalization 
condition J d 2 H,K m (R) = 1. The remaining time inte- 
gral gives the final result quoted in Eq. (|103l) . 
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